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ABSTRACT 

We present the results of a comprehensive survey using the Arecibo Observatory for Zeeman splitting 
of OH lines in OH megamasers (OHMs). A total of seventy-seven sources were observed with the 
Arecibo telescope. Of these, maser emission could not be detected for eight sources, and two sources 
were only ambiguously detected. Another twenty-seven sources were detected at low signal-to-noise 
ratios or with interference that prevented placing any useful limits on the presence of magnetic fields. 
In twenty-six sources, it was possible to place upper limits on the magnitude of magnetic fields, 
typically between 10-30 mG. For fourteen sources, the Stokes V spectra exhibit features consistent 
with Zeeman splitting. Eleven of the se fourteen are new de tections, and the remaining three are 
re-detections of Stokes V detections in iRobishaw et all (|2008[ ) . Among confident new detections, we 
derive magnetic fields associated with maser regions with magnitudes ranging from 6.1-27.6 mG. The 
distribution of magnetic field strengths suggests the magnetic fields in OH masing clouds in OHMs 
are larger than those in Galactic OH masers. The results are consistent with magnetic fields playing 
a dynamically important role in OH masing clouds in OHMs. 

Subject headings: galaxies: magnetic fields — ISM: magnetic fields — magnetic fields — masers — 
polarization 



1. INTRODUCTION 



This work follows that of lRobishaw et all (|2008f ) (here- 
after R08), who detected Zeeman splitting in a sample 
of five OH megamasers, and is structured similarly for 
ease of reference to that work. OH megamasers (OHMs) 
are powerful extragalactic masers that operate in the 
18-cm transitions of the hydroxyl molecule, with typi- 
cal luminosities of order log(LoH/^o) = 3 — one hun- 
dred million (10 s ) times greater than the luminosity of 
OH masers in the Milky Way. OHMs are found in lu- 
minous and ultra-luminous infrared galaxies (LIRGs and 
ULIRGs), classes of galaxies with log(ZfFiR,/£©) > 11 
and log(LFiR/L©) > 12, respectively. Among these 
galaxies, OHMs have a preference for the m ost far in- 
frared luminous (|Darling fc Giov anclli 2002al). 

OHM emission occurs primarily in the 1667 MHz main 
line transition, w ith the 1665 MHz main li ne either 
weaker or absent (Darling & Giovanclli 2002b), and the 
satellite lines at 1612 MHz a nd/or 1720 MHz only de - 
tected in a few near by OHMs IBaan et all (|1987l I1992D : 
iMartin et all ()1989D . VLBI observations of Arp 220, 
IRAS F17208-0014, and III Zw 35, all with redshifts 
z < 0.05, find a combination of compact 1667 MHz emis- 
sion on parsec sized scales and diffuse 1665 MHz and 
1667 MHz emission on scales of hundreds of parsecs, 
and also see evidence for Keplerian motion in rotating 
maser rings (iLqnsdale et alj|1998t [Diamond et al.H l999: 
iPihlstrom et all 120011: iMomiian et al.l I2006D . Markar- 
ian 231, meanwhile, co ntained no compac t emiss ion in 
VLBI observations by ILonsdale fc Smith! ()2003f ) that 
had a physical resolution of roughly 80 pc, meaning 
emission is coming f rom regions larger than 80 pc. 
Pihlst rdm et al.l (|2005[ ) performed VLBI observations of 
IRAS F12032+1707 and IRAS F14070+0525, which are 
the third most distant and most distant OHMs, respec- 



tively, with z > 0.20. They found that the majority of 
emission from IRAS F12032+1707 is ordered and com- 
pact on scales of less than 100 pc, while most of the emis- 
sion in IRAS F14070+0525 is resolved out. iParra et ail 
|2005T) modeled the emission of III Zw 35 as arising 
from a ro t ating starburst ring with clumpy OH clouds, 
which Lo ;2( )()•">• suggests could explain other sources as 
well. In the context o f the clumpy model of OHMs, 
lLockett fc Elitzurl (|2008l) explained the observed line ra- 
tios of OHMs as being produced through 53 /im radiative 
pumping produced by dust with a minimum temperature 
of 45 K, and overlap of lines with ~ 20 km s _1 lincwidths. 

[U] LIRGs that host OHMs are predominantly merg- 
ing systems, in which molecular gas is funnclcd to 
a central merging nucleus. One or both of star 
formation and active galactic nucleus (AGN) activ- 
ity produce far infrared and radio continuum emission 
(|Darling fc Giovanellil l2002al [III)— evidence for AGN 
or star for mation activi t y dom inating in OHMs has been 
conflicting. IBaan et al.l ()1998l ) used optical spectroscopy 
to look at the source of nuclear activity in OHMs, and 
found a preference for AGN among their sample of 42 
galaxies. IDarling fc Giovanellil ([2006) also studied opti- 
cal spectra of OHMs, this time comparing to a sample 
of non-masing [U] LIRGs, and found no significant differ- 
ence between the O HM hosts and non-masing galaxies. 
iVignali et all (|2005h observed seven OHMs with Chan- 
dra, detecting only one source weakly, a finding that is 
consistent with galax ies powered by s t ar fo rmation or 
low-luminosity AGN. Kandalian ct al. (2007) looked at 
a larger sample of x-ray observations of 22 OHMs, using 
both Chandra and XMM-Newton data, and observed a 
weak relationship between the x-ray luminosity and OH 
luminosity that favored an AGN contribution in some 
sources, but noted that galaxies may host OHMs as a 
result of either AGN activity or rapid star formation. 
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iWillett et all (|2011allbh compared OHM hosts and non- 
masing ULIRGs in the mid-infrared, finding that the 
fraction of AGN in OHMs is lower than that of non- 
masing ULIRGs, and also demonstrated support for a 
minimum dust t emperature compa r able t o that used in 
the modeling in lLockett fc Elitzurl <|2008f ) . Altogether, 
the most recent work points to star formation power- 
ing OHMs, though AGN activity may still play an im- 
portant role. In trying to distinguish O HM host s from 
non-masing galaxies of similar luminosity, iDarlingl (|2007l ) 
examined the conditions of the molecular gas, finding a 
clearer distinction between the two populations based on 
high mean molecular gas densities. Typical OHM hosts 
have n (H2) = 10 3 -10 4 cm -3 and high dense gas frac- 
tions, with half of galaxies that have ^hcn/^co > 0.07 
also hosting OHMs. 

R08 noted that the conditions in which OHMs reside, 
particularly the high energy and molecular gas densities, 
suggest the presence of strong magnetic fields. Minimum 
energy arguments made using VLA observations of syn- 
chrotron emission argue for minimum volume averaged 
magnetic field strengths of 1 m G in the central 100 pc o f 
ULIRGs (jCondon et al.1 119911 ). iThompson et all (2006) 
argued that producing the FIR-radio correlation for star- 
burst galaxies requires significant fine-tuning if magnetic 
fields are not significantly larger than the minimum en- 
ergy estimates. They further show that for magnetic 
fields to be dynamically important relative to gravity, 
as in the Milky Way, the magnetic fields in Arp 220 or 
a similar galaxy would be of order 30 mG. These argu- 
ments for fields in [U]LIRGs with magnitudes between 1- 
30 mG provided a compelling case to search for Zeeman 
splitting in [U]LIRGs, and in particular, observe OHM 
lines. The hydroxyl molecule has a relatively large mag- 
netic dipole moment, and OHMs provide bright, narrow 
lines in which Zeeman splitting may be detected. 

In an effort to directly Zeeman splitting and measure 
line-of-sight magnetic fields in [U]LIRGs, R08 used the 
Arecibo telescope and the Green Bank Telescope to per- 
form Full-Stokes observations for eight OHMs. They de- 
tected magnetic fields with a median strength of ~ 3 mG 
in five of the eight observed galaxies, and the strongest 
detected field had a magnitude of 17.9 mG. Their re- 
sults confirmed that Full-Stokes observations of OHMs 
provide a viable way to directly measure magnetic fields 
in [U]LIRGs, even if magnetic fields are not of equal dy- 
namical importance to gravity in the ISM of [U]LIRGs. 
Here, we go further, observing seventy-seven known ex- 
tragalactic OH masersQ including re-observations of six 
of the sources from R08. Eight OHMs are not detected at 
all, primarily as a result of increased radio-frequency in- 
terference (RFI) since their discovery. Three sources ex- 
hibit features in their Stokes V spectra consistent with 
Zeeman splitting, but in which we are only marginally 
confident. Eleven OHMs have Stokes V features that we 
confidently associate with Zeeman splitting, including re- 
detections of Zeeman splitting in three OHMs observed 
in R08. One source, IRAS F10173+0829, showed ap- 
parent Zeeman splitting in the observations of R08, but 
re-observation showed that the apparent features were 

1 While most meet the definition of megamaser, a few would 
more accurately be called kilomasers or gigamasers. For simplicity, 
we subsequently call all of these sources OHMs. 



interference. The remaining sources were detected in 
Stokes /, but did not show significant Zeeman splitting 
detections. 

2. SOURCES 

The seventy-seven OHMs make up the entire sam- 
ple known at the time of observations that are ob- 
servable with the Arecibo telescop43 in Puerto Rico. 
All but two of the targets were taken from the OHM 
su rvey and compilation of p a st de t ections presented 
in iDarling fc Giovanellil (|2000l [200ll l2002a[) (hereafter 
DG00, DG01, DG02). The two exceptions are the only 
two Arecibo accessib le sources tha t have been discovered 
in the decade since (|Willettll2012t l. While the choice to 
use only Arecibo excludes some viable targets, the re- 
sults of R08 showed that detecting magnetic fields using 
Zeeman splitting of OHMs requires a level of sensitivity 
more easily achieved with Arecibo. 

3. OBSERVATIONS 

We used the L-band wide receiver on the 305- 
meter Arecibo telescope in full-Stokes mode to observe 
all known OH megamasers at declinations accessible 
to Arecibo, over a period spanning December 2007- 
December 2009. The interim correlator was set up to 
observe four different bands: one 12.5 MHz band, cen- 
tered halfway between the 1665 and 1667 MHz, and three 
6.25 MHz bands, centered at 1612, 1667, and 1720 MHz. 
This set-up allowed simultaneous observation of all four 
ground state OH maser lines and all four Stokes param- 
eters, /, Q, U, and V. This differs from the correlator 
setup used by R08, which focused only on the main lines, 
and had a range of bandwidths that depended upon the 
velocity extent of the source. The uniformity of the setup 
used here eased reduction and analysis. This work fo- 
cuses on the 1665 MHz and 1667 MHz lines, as there were 
no Zeeman detections in the satellite lines at 1612 MHz 
or 1720 MHz. A future work will present the results of 
the observations of the satellite lines. 

We spent equal time on and off source, switching po- 
sition every 4 minutes. The off source position had the 
same declination and a right ascension 4 minutes east 
of the source, which kept the hour angles of on and off 
source observations nearly equal. The integration time 
was 1 second, chosen to allow for elimination of RFI that 
appeared on short time scales. 

4. DATA REDUCTION 

Many of the details of the data reduction are well de- 
scribed in R08. This includes calibration of spectra, RFI 
removal, bandpass and gain correction, and fitting Gaus- 
sian components to the Stokes I profiles. The discussion 
here will focus on describing new components of the data 
reduction, and mostly does not repeat the description in 
R08 in cases where the procedure used in reducing the 
data in this paper was identical to that in R08. Section 
4.1 of R08, describing calibration, and 4.3, describing 
bandpass and gain correction, are omitted entirely. 

2 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 
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We use the classical definition of Stokes /, which is the 
sum, rather than the average, of two orthogonal polar- 
izations. Our flux densities are thus twice those of most 
previously published work, but are consistent with the 
choice of R08. We take the IAU definition of Stokes V, 
with V = RHCP - LHCP. Following the IEEE defini- 
tions RHCP is clockwise rotation and LHCP is coun- 
terclockwise rotation of the electric field vector along the 
direction of propagation. 

4.1. RF I removal 

The data were examined for RFI in various stages. 
First, the on/off pairs of 1 second records were flagged 
for RFI and checked for significant rippling caused by 
the Sun or Moon being within the beam sidelobes. At 
this point, the data reduction forked: within each on/off 
pair of 4 minute records, the 1 second records were av- 
eraged together to form a single 4 minute spectrum, and 
also combined by taking the median to form a second 
4 minute spectrum for the same data. In cases where the 
RFI is limited, the average spectrum tended to be less 
noisy, but forming a composite spectrum from the me- 
dian of the 1 second integrations was clearly better for 
sources with significant narrow-band RFI. 

Each source was typically observed for a few hours. 
This resulted in ^10-20 4 minute on and off spectra for 
each source. These too were combined in the same man- 
ner as the 1 second spectra, either by taking the average 
or the median. Applying these choices to each set of 
mean or median 4 minute spectra resulted in four differ- 
ent final spectra for each source at a given frequency and 
polarization: mean/mean, mean/median, median/mean, 
median/median. Of these four options, the spectrum 
with the least noise was selected. 

4.2. Confidence tests 

Many of the Zccman measurements presented here 
represent marginal detections. We describe here some 
of the primary obstacles to determining magnetic field 
strengths, and tests applied to the data and fits that val- 
idate the results. 

4.2.1. Broadband rippling 

Many of the Stokes V spectra contain ripples of char- 
acteristic size ~ 0.3-1 MHz that can masquerade as Zee- 
man splitting of broad OHM components, and represent 
the most serious source of systematic error in estimat- 
ing magnetic fields. The best fit fields to these ripples 
in the Stokes V spectra are often in excess of 50 milli- 
gauss. These ripples are easily spotted by eye in most 
situations, but can confuse magnetic field fits in sources 
with particularly rich spectra. For this reason, magnetic 
field fits associated with spectral features wider than 0.3 
MHz are all considered to be untrustworthy. This choice 
is also physically motivated, in that particularly wide 
maser lines are likely produced by blending from many 
masing clouds, and are less likely to provide a detectable 
line-of-sight Zeeman splitting signal. 

4.2.2. Spectral shifting 

Even after removing the broadest ripples, there are 
Stokes V spectra with structured noise on smaller fre- 
quency scales that can confuse interpretation. In some 



cases, noisy structure in Stokes V may be coincidentally 
lined up with features in the Stokes I spectrum, and 
produce claimed magnetic field detections. To identify 
instances of dubious magnetic field claims, we ran our 
fitting procedure for Zeeman splitting with the Stokes V 
spectrum shifted relative to the Stokes / spectrum for all 
sources. With 2048 channels, this produces 2048 differ- 
ent fits to features in Stokes V. For each set of fits, we 
selected the features narrower than 0.3 MHz, found the 
associated signal to noise ratio for the magnetic fields of 
each narrow component, and produced an overall quality 
of fit at each shift. For sources where non-zero shifts pro- 
duced fits of comparable quality to the fits with no shift, 
the reported fields were regarded as highly dubious. In 
Section [5X71 all sources for which we claim a magnetic 
field detection were best fit when no shift was applied 
between the Stokes / and V spectra. Other sources with 
claimed magnetic fields, but which were also well fit when 
the Stokes I and V were shifted relative to one another, 
were considered Stokes V non-detections or marginal de- 
tections, depending on other factors. An example of this 
process is shown in Figure [T] The source in the top panel 
is not fit particularly well at any one shift, whereas the 
source in the bottom panel is only well fit when no shift 
is applied, providing confidence that the features being 
fit are real. The full fits for these sources are shown in 
Sections 15.1.51 and 15.1.71 

4.2.3. Comparison of overlapping spectra 

For each source observed, there was a total of 25 MHz 
of spectral coverage, with two 6.25 MHz bands cover- 
ing the OH satellite lines at rest frequencies of 1612 and 
1720 MHz, and the broad 12.5 MHz band that encom- 
passed both main lines at 1665 and 1667 MHz. With 
sources spread out over a redshift range of 0.007-0.265, 
the highest observed frequency was 1711 MHz, and the 
minimum was 1271 MHz. Just over 90% of this range 
was covered by at least one source, including the entire 
range between 1351.8 MHz-1663.7 MHz. Moreover, 73% 
of this band was covered by at least two sources, and 
50% was covered by at least four sources. 

This overlapping coverage was a valuable tool in RFI 
identification. Spectra with features of unclear origin 
could be compared to spectra of other sources in the 
same frequency range. If a similar feature appeared 
in both spectra, then it strongly suggests interference. 
Comparing overlapping spectra of different sources was 
critical in discovering that the Stokes V features of 
IRAS F10173+0829, reported in R08, were actually in- 
terference. A rough characterization of RFI severity 
across the frequency range in the survey can be seen in 
Figure [2 

4.2.4. Field sign 

In a survey of Zee man splitting in Galactic OH masers, 
Fish & Reid (2006) tested for unknown systematics in 
their data by counting how many Zeeman pairs had 
greater flux in the LCP component, how many had 
greater flux in the RCP component, and comparing their 
result to the expectation that the each possibility was 
equally likely. We performed a similar test, counting the 
number of positive and negative magnetic fields derived, 
and comparing the results with the binomial distribution. 
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Fig. 1. — The average signal-to-noise ratios of fitted fields 
to narrow components in the Stokes / and V spectra of IRAS 
F23028+0725 (top panel) and IRAS F09039+0503 (bottom panel) 
are plotted for fits when the Stokes I and V spectra are shifted 
relative to each other. The overall fit quality at no shift for 
F23028+0725 is mediocre, and the fits are marginally better when 
the spectra are shifted relative to each other, suggesting that the 
features being fit in the Stokes V spectrum are not Zeeman split- 
ting. The fits for F09039+0503 are at their best for no shift, and 
no other shift produces fits of comparable quality. 



For N fields and an equal likelihood of positive or neg- 
ative fields, the probability of a positive field is p = 0.5, 
the expected number of positive fields is Np, the ex- 
pected number of negative fields is N(l — p), and their 
variance is Np(l —p). For all sources, we fit a total of 373 
Gaussian components, and a magnetic field estimate was 
made for each component, giving an expected number 
of f86.5 ± 9.7 for positive and negative fields. Of de- 
rived fields, 172 had positive sign, and 201 were negative 
in sign, which are not statistically significant deviations. 
Likewise, for magnetic field components we considered a 
detection, there were a total of 35 derived fields. Com- 
pared to an expectation of 17.5 ± 2.1, 19 were positive 
in sign and 16 were negative, very close to the expected 
distribution. 

4.2.5. Bootstrap resampling 

An often used technique to distinguish a source of 
broadband interference from one of astrophysical interest 
is to split a dataset in half, and check that the feature is 
roughly the same strength throughout the observation. 
When looking for very low SNR features, this method 
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Fig. 2. — The spectral coverage of the survey is shown here, along 
with the RFI environment at Arecibo during observations. Areas 
with observations and no observed RFI are plotted at a default 
level of 1 mjy, as that was the typical rms error for RFI free data. 
At frequencies where RFI was present in at least one data set, 
the approximate amplitude of the most severe RFI is shown. For 
a couple of wide regions with severe RFI, such as that between 
1524—1552 MHz, the value plotted is typical of the large RFI spikes 
throughout the frequency range, but RFI of that amplitude was not 
necessarily observed in each channel. 



proved to be only marginally helpful, as the fluctuations 
in amplitude for a low SNR astrophysical signal due to 
noise might be similar to those expected from broad- 
band interference. As an extension of this, we bootstrap 
resampled our spectra, rather than splitting the data in 
to two halves. 

Bootstrapping is a resampling method that can be 
used for estimating confidence inter vals associated with a 
statistic (|Efron fc T ibshirani 1993)). For a statistic such 
as the median, bootstrap resampling makes it possible 
to characterize uncertainty, despite there being no formal 
way to do so. The procedure is relatively straightforward. 
Given a sample of n observations, x = (x\, X2, x n ), the 
empirical probability of observing each of the n values is 
1/n. A bootstrap sample of x is made up of n values 
drawn randomly from x. One of the original observa- 
tions may be drawn more than once, or not at all, in 
this bootstrap sample. To build a distribution of values 
for some statistic, such as the median, take the median 
of the original sample as well as the median of each of 
many bootstrap samples. This sampled distribution of 
the median measurements provides an estimate of the 
uncertainty in the measurement of the median. 

As the spectrum for each OHM was made by combin- 
ing 4 minute spectra, either using the mean or median, 
we treated each of the 4 minute spectra as an individ- 
ual observation. As discussed above, for some sources, 
taking the median of those 4 minute spectra produced a 
spectrum with less noise, usually because it better han- 
dles RFI, while for some sources the mean spectra were 
less noisy; the "measurement" in this case is whichever of 
the two was better. We then took 1000 bootstrap sam- 
ples of that measurement. The resultant lower and upper 
limits are in a few cases useful in establishing whether a 
broad spectral feature is OH emission. We performed 
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similar resampling in checking the derived errors of mag- 
netic held estimates. For sources with claimed mag- 
netic field detections, the Stokes I and Stokes V spectra 
were jointly resampled, and the magnetic field deriva- 
tion procedure described below was applied. The error 
from bootstrap resampling was consistently comparable 
to those reported in the fitting procedure. 

4.3. Fitting Gaussian Components to Stokes I Profiles 

We followed the same general guidelines given in R08 
for fitting Gaussians to the overall Stokes / profile. Fits 
were made by eye, so are inherently subjective, but fol- 
low relatively straightforward guidelines that aim to pro- 
vide a balance between minimizing the number of Gaus- 
sian components while also accounting for the majority 
of structure in the OHM profiles. 

5. RESULTS 

5.1. Circular Polarization and Line-of-Sight Magnetic 

Fields 

5.1.1. Magnetic Field Derivation 

As in the section on data reduction, much of the pro- 
cess for obtaining magnetic field measurements is the 
same as R08. In particular, magnetic fields for most 
sources are derived assuming that the Zeeman splitting 
produced by the magnetic field is smaller than the intrin- 
sic linewidth. In this situation, the Stokes V spectrum is 
proportional to the derivative of the Stokes I spectrum, 
with 



V 



cl. 



(1) 



where Bu is the line-of-sight magnetic field, b is the split- 
ting coefficient (jHeiles et al.l Il993h . and the cl term ac- 
counts for leakage of / into V. The splitting coefficient 
is itself proportional to the Lande g-factor for the ob- 
served transition; for the 1667 MHz line in OH, it is 
1.96 Hz fiG- 1 , and for the 1665 MHz OH line, it is 
3.27 Hz /iG -1 . In deriving line-of-sight magnetic fields, 
we perform least squares fits to Equation!]] As explained 
and justified in R08, the fits to the Stokes I profile serve 
as the independent variable in solving Equation [TJ For 
each Gaussian component used in the overall Stokes / 
profile, a separate line-of-sight magnetic field is fit, along 
with estimated errors. The reported errors are purely 
statistical, and underestimate uncertainty in many spec- 
tra. For that reason, we do not rely on the reported 
errors alone when assigning confidence to a fitted mag- 
netic field. 

This approach does not account for the blending of 
magnetic field measurements that may occur if multiple 
masing clouds contribute to the observed emission at a 
specific velocity, which likely intr oduces error for at l east 
some fraction of observed OHMs. iParra et al.l ()2005[ ) ex- 
plained the combined compact and d i ffuse emission in 
III Zw 35 observed by Pihlstrom etH] ([2001 as a result 
of clumpy OH masing clouds in a narrow ring. Bright, 
compact maser features are produced when more than 
one cloud is along the line-of-sight. R08 compared their 
single dish spectru m of III Zw 35 w i th the interferomet- 
ric observations of IPihlstrom et all (|2001l) . They found 
that three of the Gaussian components that they fit to 



the total III Zw 35 profile correspond to t hree compact 
features mapped by IPihlstrom et ail (|200lD . 

This adds a layer of complication to interpretation of 
Zeeman splitting results, as well as making detections 
more difficult. For instance, two masing clouds at the 
same velocity along the line-of-sight that have equal in- 
tensities and magnetic fields of equal magnitude but op- 
posite sign will combine to produce a perfectly flat Stokes 
V spectrum, meaning no magnetic field would be mea- 
sured. More generally, the measured magnetic field is the 
average of the intensity weighted fields for each cloud. 
Thus the reported fields may reflect lower magnetic field 
magnitudes than are actually being probed, if blended 
lines are associated with fields oriented in opposite di- 
rections. 

Descriptions of each source follow. Described prop- 
erties of OH megamaser host galaxies are bas e d upo n 
the results compiled in IDarling fe Giovanellil 12 000). 
iDarling fe Giovanellil (|2001|L and IDarling fe Giovanellil 
(|2002al) . Lines are described as "narrow" or "wide" 
throughout the following sections. Wide lines are those 
fit with Gaussians that have a full width at half maxi- 
mum (FWHM) greater than 0.3 MHz, and narrow lines 
have FWHM less than 0.3 MHz. This cutoff is motivated 
by the size of noisy rippling in Stokes V spectra that can 
masquerade as Zeeman splitting, making reported fields 
on wide lines highly dubious. For sources with features 
consistent with Zeeman splitting, and a small sample of 
sources with no Zeeman splitting, we present a table with 
the full fitting results. Narrow lines with fitted fields that 
are > 3er are highlighted in bold text. 

5.1.2. Re-observed sources 

The sources with Full-Stokes observations made by 
R08 in February 2006 were observed again in 2008-2009 
as part of this survey. This was done for two reasons: 
to look for evidence of variability in OHMs, and to con- 
firm the magnetic field strengths derived by R08. Over- 
all, our results are very similar to those presented in 
R08, though our more recent observations are generally 
slightly noisier, and have higher reported errors associ- 
ated with Gaussian components and fields. For this rea- 
son, we only reproduce the table with fits for one source, 
IRAS F12032+1707, in order to discuss VLBI observa- 
tions. However, we note that for all sources presented in 
R08, the tables reported an intensity for Gaussian 4 that 
is too small by a factor of 100. 

We find evidence of weak variability in the Stokes / 
profiles of most of the sources, and no strong evidence of 
variability in the Stokes V profiles. For one source, IRAS 
F10173+0829, we discovered that apparent Zeeman split- 
ting seen in R08 was actually circularly polarized RFI. 

IRAS F01417+1651 (III Zw 35): The shape of fea- 
tures in the Stokes I spectrum remain si milar in the more 
recent observations, as shown in Figure 3(a) but the en- 
tire 1667 MHz feature appears to be moderately brighter 
in 2008-2009 than it was in 2006. The peak flux density 
feature, located at 1622.8 MHz, brightened by ^40 mJy, 
a fractional change of ~8%. There are no clear changes 
between the old and more recent Stokes V spectra. The 
individual Gaussian components used for fitting the new 
spectrum, and the magnetic field values associated with 
each, are in good agreement with those of R08. In addi- 
tion, the same field reversal is observed. Despite bright- 
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(c) IRAS 12112+0305: The difference between the two spectra (d) IRAS F15327+2340: The difference between the Stokes J 
is offset from a mean close to zero. profiles was scaled by a factor of 4. 



Fig. 3. — Stokes / (top panel) and V (bottom) profiles four sources observed first in 2006 and then again in 2008-2009. The data taken 
in 2006 was first presented in R08, and shown here again with dashed lines, while the more recent data are shown with a thin solid line. 
The difference between the two data sets (thick solid) are in some cases offset or scaled to show variation more clearly. 
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ening, the structure of the OHM in IRAS F01417+1651 
appears relatively unchanged. 

IRAS F10173+0829: R08 reported features consis- 
tent with Zeeman splitting in their observations of this 
OHM. Unfortunately, the Stokes V spectrum composed 
by taking the median of 1-second and 4-minute spec- 
tra found weaker features in the Stokes V spectra than 
were seen by R08. Further, the Stokes V spectra of other 
sources at the same frequency had Stokes V features sim- 
ilar to those seen in IRAS F10173+0829. Figure H shows 
the Stokes V spectrum of the 1612 MHz line of IRAS 
F15107+0724 compared to that of IRAS F10173+0829. 
IRAS F15107+0724 showed no evidence of maser emis- 
sion at frequencies between 1589.2-1589.4 MHz in its 
Stokes I spectrum, yet had a similar Stokes V signal to 
IRAS F10173+0829'. This not only means that the ap- 
parently detected Zeeman splitting is not real, but also 
prevents placing meaningful limits on magnetic fields in 
this OHM. 
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Fig. 4. — The solid line shows the Stokes V spectrum for 
IRAS 10173+0829, while the dashed line shows the Stokes V 
spectrum for IRAS 15107+0724, which is of similar amplitude 
and shape, at the same frequency. While the features in 
IRAS 10173+0829 nicely line up with maser emission in its Stokes 
I spectrum, the Stokes / spectrum for IRAS 15107+0724 does not 
have any features suggesting maser emission, indicating that the 
Stokes V features in both sources are interference. 



IRAS F12032+1707: The results of the more recent 
observations are in gener al agreement with those of R08, 
as shown in Figure 3(b) The rms error in the difference 
spectrum is dominated by that in the more recent spec- 
trum, which has ~2 mJy channel-to-channel variation. 



There is weak evidence for variation associated with the 
brightest, narrowest features in the Stokes / spectrum, as 
the difference spectrum has small peaks with amplitudes 
of ^5 mJy associated with each bright emission feature, 
but this variation is nowhere near as dramatic as the 
nearly factor of two brightening observed in the brightest 
peak, located at 1371.4 MHz, between the observations 
of DG01 and those of R08. The Stokes V features remain 
consistent between observations, and the magnetic fields 
found by fits to the new spectra are entirely consistent 
with those published in R08. 

There are exis t ing V LBI observations of this OHM in 
IPihlstrom et al.l ((20 05) that were not discussed in R08. 



For this reason, Table [T] is reproduced from that paper. 
The emission is compact on scales of less than 100 pc, and 
all of the single dish flux is recovered in this region, within 
error. There is a clear north-south veloci ty gradient, with 
the bluest peak located in the north. IPihlstrom et al.1 
(2005) noted a peak centered at 64,723 kms -1 that 
appeared asymmetric in the R08 observations. R08 fit 
this peak with two Gaussian components, at velocities of 
64,737 km s -1 and 64,715 kms -1 . The Gaussian com- 
ponents, when used to fit the Stokes V spectrum, yielded 
fields of 10.9±1.7 mG and 17.9±0.9 mG, respectively. 

IRAS F12112+0305: In comparing t he ol d and new 
spectra for this OHM, shown in Figure 3(c) the most 
obvious difference is the appearance of weak interference 
spikes in the older spectra that were not a problem in 
the more recent observations. There is also astrophysi- 
cal change, as the brightest feature in the spectrum for 
IRAS F12112+0305 has apparently dimmed by ~10 mJy, 
relative to a flux density of ~100 mJy. At the same 
frequency as the change in the Stokes / spectrum, the 
Stokes V spectrum has also changed, as there is now a 
broad dip associated with the peak of the maser emission. 

Despite this new dip, combined fits to the Stokes / 
and V profiles do not produce any evidence of magnetic 
fields. At the frequencies corresponding to 1667 MHz 
emission, the more recent spectra have comparable noise, 
and errors in derived Gaussian fits to the structure, as 
the original observations. We find a comparable upper 
limit of ^3 mG for magnetic fields in the masing regions. 

IRAS F14070+0525: R08 were only able to only 
place upper limits on magnetic fields in this OHM. Our 
Stokes I and V spectra are roughly a factor of two noisier 
than those presented in R08. Consequently, no Stokes V 
features are detected, and the upper limits we may place 
on magnetic fields are less stringent. Comparing the old 
and new spectrum, there is no evidence for variations 
with amplitudes greater than 3 mJy. 

IRAS F15327+2340 (Arp 220): Vari ability has pre- 
viousl y been obs erved in Arp 220 by lLonsdale et al.l 
(2008). In Figure 3(d)| our spectra are compared with 
those of R08, and show strong support for claims of vari- 
ability. In the Stokes / spectrum, narrow features be- 
tween 1637.5-1638.3 MHz vary by 10-15 mJy between 
the old and new spectra, some dimming and some bright- 
ening, as well as a slight brightening of the "shoulder" 
at 1638.55 MHz. Given the brightness of Arp 220, these 
difference represent only ~1% variations in the flux den- 
sity. To the eye, the most noticeable difference in the 
Stokes V spectrum is a ~0.5 MHz hump that falls in 
the same frequency range as the changes in the Stokes I 
structure, but it is difficult to make out clear differences 
between the narrower peaks. 

Overall, when applying the same fitting routine to the 
newer spectra and comparing to the older spectra, the 
features are largely the same, with fitted components 
having center frequencies, widths, amplitudes and associ- 
ated magnetic fields that are the same within error of the 
fits. The one difference that does appear is in a Gaussian 
at 1638.3 MHz, which in the original spectrum fit with a 
magnetic field with strength 2.03 ±0.76 mG, whereas the 
more recent spectrum was fit with a field with strength 
7.29 ± 1.43 mG. Upon careful inspection, the more re- 
cent Stokes V spectrum is deeper at that frequency, and 



TABLE 1 

IRAS F12032+1707 Gaussian Fit Parameters 





S 


V 


Au 




Vq 




Gaussian 


(mjy) 


(MHz) 


(MHz) 


(km s~ - 1 ) 


(mG) 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 





4.31 ± 0.54 


1367.1002 ± 0.0169 


0.4082 ± 


0370 


65844.0 


— 6.79 ± 7.94 


1 


3.95 ± 0.96 


1367.6739 ± 0.0457 


0.7915 ±0 


1793 


65690.6 


30.42 ± 12.43 


2 


8.42 ± 0.99 


1368.8016 ± 0.0094 


0.6264 ±0 


0448 


65389.5 


6.91 ±4.86 


3 


3.00 ±0.34 


1369.5922 ± 0.0071 


0.1460 ±0 


0207 


65178.7 


-1.05 ±6.48 


4 


0.21 ±0.36 


1369.6903 ± 0.0452 


2.6434 ± 


1574 


65152.6 


- 3.25 ± 3.98 


5 


3.23 ±0.54 


1369.9216 ± 0.0036 


0.0450 ± 


0089 


65091.0 


3.10 ±3.41 


6 


6.67 ±0.34 


1370.7194 ± 0.0064 


0.2954 ±0 


0165 


64878.6 


11.64 ±4.31 


7 


21.67 ±0.65 


1371.2516 ± 0.0030 


0.4769 ± 


0115 


64737.0 


10.90 ± 1.72 a 


8 


16.99 ±0.46 


1371.3316 ± 0.0010 


0.0849 ± 


0029 


64715.8 


17.92 ±0.89~ 


9 


16.19 ±3.38 


1372.1576 ± 0.0568 


0.4872 ±0 


0531 


64496.4 


1.78 ± 2.48 


10 


25.02 ±4.51 


1372.3299 ± 0.0027 


0.2795 ±0 


0165 


64450.6 


- 1.45 ± 1.12 


11 


9.56 ±0.51 


1372.7780 ± 0.0412 


0.7315 ±0 


0588 


64331.7 


-11.69 ±4.98 


12 


2.02 ±0.21 


1373.8371 ± 0.0179 


0.3807 ±0 


0500 


64051.0 


20.61 ± 15.49 



a Derived magnetic field strengths above 3cr that are considered believable are marked in bold. 



the difference has an amplitude larger than the features 
in the difference spectrum at frequencies not associated 
with the Stokes / spectrum. Nonetheless, given other 
amplitude variations in the Stokes V spectrum nearby, 
we do not consider this strong evidence of variation in 
the magnetic field strength. 

5.1.3. Stokes I non- detections 

Eight sources previously reported as OHM hosts were 
not detected in our survey. In most cases, non-detections 
are due to increased RFI since the original discovery of 
these sources. 

For non-detected sources, upper limits are calculated 
according to the prescription in DG02, with 

This assumes a boxcar line profile of height 1.5 a, which 
is the RMS noise in the Stokes / spectrum, and a rest 
frame velocity line with of 150 km s _ , which is the av- 
erage FWHM of the 1667 MHz line of the known OHM 
sample. To calculate the luminosity distance, Z?l, we use 
the results from WMAP5 and assume Ho = 70.5 km s -1 
Mpc" 1 , n A = 0.726, and ti M = 0.274 (jHinshaw et al.l 
2009). The factor of i appears because we use the clas- 
sical definition of Stokes I, equal to the sum of the or- 
thogonal polarizations. 

In some cases, comparison with previo u sly published 
work is diffic ult. IBottinelli et"aH (il987l) . IMartin et all 
(|1989fh and IBottinelli et al.l (|1990f) "do not explicitly 
state the cosmology used to calculate the isotropic lumi- 
nosities they reported. For comparison wi th the values 
reported in those instances, we follow IMartin et al.1 
(1198811 . who used H Q = 75 km s _1 Mpc -1 and a 
deceleration parameter, qq = 0. 

I RAS F00509-I - 1225: This OHM was first reported 
by IBottinelli et~aT1 (jl990h . with log(L OH /i©) = 1-78, 
but has no published spectrum. Given the reported red- 
shift Zq — 0.060875, the li kely luminosity dista nce they 
used is 243 Mpc. For the (jHinshaw et all 12009ft cosmol- 
ogy, the luminosity distance at this redshift is 271 Mpc, 
implying an isotropic luminosity of log(ioH/^o) = 1.88. 



At this redshift, maser emission should be expected near 
1571 MHz. While the composite spectrum has features 
near this frequency, those features completely disappear 
or even appear as absorption in bootstrap resampling, 
suggesting they are instead due to interference from GPS 
L2 that is centered at 1575.42 MHz. The rms error is 
5 mjy, so we place an upper limit log(LoH/^o) < 1-82 
on the isotropic luminosity of OHM emission, weakly sug- 
gesting that this source is n ot actually an O H M. 

IRAS F03056+2034: IBottinelli et~all (|1990f ) re- 
ported this OHM with a an isotropic luminosity of 
log(£oH/£©) = 1-28. Located at a redshift of z — 
0.027369, the expected maser emission is in a frequency 
range dominated by spiky interference likely from the 
Iridium satellite, so much so that even the median com- 
posite spectrum had an rms of 15 mjy. From this, we 
place an upper limit on the isotropic luminosity of the 
OHM emission that we could detect of log(ioH/^o) < 
1.60. This is larger than the previously reported value, 
so the non-detection is unsurprising. 

IRAS F11069+2711: No spectrum of this OHM is 
available, nor are there an y published p roper ties. It was 
listed as a megamaser by iBaan et al.1 (fl9 98?) . There is 
spiky interference in our spectra, so we took the me- 
dian of sho r ter sp ectra to produce the total spectrum. 
IBaan et all (|1998ft assumed a redshift for this source of 
of z — 0.0703, which places the expected frequency for 
OHM emission at 1557.8 MHz. There is a roughly 5 mjy 
feature centered around 1557.5 MHz, but this is likely 
due to terrestrial interference, as it appears intermit- 
tently in spectra both on and off sourc e. A more re- 
cent redshift measure by lLawrence et al.l (|1999f) puts the 
galaxy at a redshift of z = 0.072971, shifting the expected 
frequency for OHM emission to 1554 MHz. At this fre- 
quency, there is nothing resembling emission. Outside of 
the interference at 1557.5 MHz, there is an rms error of 
2.3 mjy, which we use to place an upper limit on the 
OHM luminosity of log(L H/i©) < 1-65. 

IRAS F13451+1232: With a 5.4 Jy continuum flux 
at 1.4 GHz, this OHM is the strongest continuum source 
in our sample. A te ntative OH detection was claimed by 
iDickev et al.l (Jl990), who reported a flux of 1.7 mjy and 
an isotropic luminosity of log(LoH/£©) = 2.38. While 
most spectra in our sample have rms errors of order 
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(a) IRAS 08071+0509: Zeeman splitting non-detection. 
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(b) IRAS F10339+1548: Zeeman splitting non-detection. 



Fig. 5. — Each sub-figure shows Stokes I and V spectra in two panels. Sub-figures (a), (b), and (c) show non-detections of Zeeman 
splitting. Sub-figures (d), (e), and (f) are considered marginal Zeeman splitting detections. The remaining sub-figures are new Stokes 
V detections. In each sub-figure, the top panel shows the Stokes I spectrum for each source. The flux density is defined as the sum of 
the orthogonal polarization, / = RHCP + LHCP, rather than the average. The data are shown by the solid line. The dashed lines are 
Gaussian components used to fit the overall profile. The thin solid line shows the residuals of the fits, expanded by a factor of three. The 
bottom panel in each sub-figure shows the Stokes V profile. The data are indicated by the solid line. The fit, produced using the Gaussians 
in the top panel as inputs, is shown with a dashed line. 



1 mjy, the combination of strong continuum emission 
and interference in the spectrum for this source resulted 
in an rms error of ~ 20 mjy, so we can neither confirm 
nor refute this tentative d etection. 

IRAS F15233+0533: iBaan et all (|1998ft listed this 
as a megamaser, but no published spectrum or proper- 
ties are available. At a redshift of z = 0.054064, OHM 
emission is expected at 1581.8 MHz, which is near in- 
terference from GPS L2 located at 1575.42 MHz. There 
are no features in the final spectrum that are clearly per- 
sistent. Bootstrap resampling of the 4-minute spectra 
shows that features in the composite spectrum are not 
clearly persistent. With an rms error of 20 mjy, we place 
an upper limit of log(LoH/£©) < 2.33 for OHM emission 
from this galaxy. 

IRAS F17526+ 3253: This sou r ce wa s first observed 
for OH emission bv lGarwood et al.1 (|l987i), who reported 
a non-detection with a spectrum that had 1.0 mjy 
rms noise, which corresponds to an upper limit on the 
isotropic luminosity of log(LoH~/ Lq) < 0.60 for the cho- 
sen cosjnologj^Aii_OHM was reported two years later 
by iMartin et al.l (|T989), with an isotropic luminosity of 
log(LoH~/ Lq) = 0.99, but no accompanying flux or spec- 
trum was published. In our spectrum, interference from 



the Iridium satellite is severe. After combining the me- 
dian spectra, the rms error is ^10 mjy, so we do not 
detect OHM emission, and cannot place any useful new 
limit on OHM emission from this galaxy. 

I RAS F20491+1846: This OHM was first reported 
bv IBottinelli et alj (|1989[ ). There is no published flux, 
and the published isotropic luminosity is log(LoH/^o) = 
1.09 at a redshift of 0.0290. Even looking at the me- 
dian spectrum, spiky interference from the Iridium satel- 
lite dominates. With roughly 30 mjy rms error in the 
spectrum, the upper limit on the isotropic luminosity is 
log(LoH~/-k©) = 1-96, well above the reported luminosity. 

IRAS F23135+2517: Interference from the Iridium 
satellite is present in the spectrum for this source. Com- 
bining the median of the me dian spectra, the rms error 
in the spectrum is ~ 20 mjy. IMirabel fc Sandersl (fl987h 
reported a flux of 2.4 mjy for this source, corresponding 
to 4.8 mjy for our definition, so the OHM is not detected. 

5.1.4. Ambiguous Stokes I detections 

IRAS F15250+ 3608: This OHM was reported in 
an IAU Circular by IBottinelli et al.1 (|1987l) to have an 
isotropic luminosity of log(LoH/-k©) = 2.58. At a red- 
shift of z = 0.0554, OHM emission is expected at roughly 
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(c) IRAS F23028+0725: Zeeman splitting non-detection. 



(d) IRAS F15224+1033: Marginal detection of Zeeman splitting. 



Heliocentric Velocity (km s ') 
41500 41000 40500 




1464 1465 1466 1467 1468 1469 1470 
1 i 3 4 S 67 8 }' 11 12 




1464 1465 1466 1467 1468 1469 1470 
Observed Frequency (MHz) 



Heliocentric Velocity (km s ') 
11200 11100 11000 10900 10800 10700 10600 



120 
100 

1 60 
40 
20 




'I I ' 



'I I 



20550+1655 



3 mjy 




1607.0 1607.5 1608.0 1608.5 1609.0 1609.5 1610.0 1610.5 

i3',6 7 8 9 10 



E 
> 




1607.0 1607.5 1608.0 1608.5 1609.0 1609.5 1610.0 1610.5 
Observed Frequency (MHz) 



(e) IRAS F15587+1609: Marginal detection of Zeeman splitting. (f) IRAS F20550+1655: Marginal detection of Zeeman splitting. 



Fig. 5. — (continued) 
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(g) IRAS F02524+2046: Zeeman splitting detection. 
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(h) IRAS F04332+0209: Zeeman splitting detection. Instead of 
showing fits to Stokes I, the RHCP (dashed) and LHCP (dash- 
dot) are plotted in the top panel, and RHCP — LHCP for the fitted 
Gaussian is shown in the bottom panel. 
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(i) IRAS F09039+0503: Zeeman splitting detection. 
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(j) IRAS F10378+1108: Zeeman splitting detection. 



Fig. 5. — (continued) 
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(k) IRAS F16255+2801: Zeeman splitting detection. 
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(1) IRAS F18368+3549: Zeeman splitting detection. 
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(m) IRAS F18588+3517: Zeeman splitting detection. 
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(n) IRAS F22134+0043: Zeeman splitting detection. 



Fig. 5. — (continued) 
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1580 MHZ, near interference from GPS L2 at 1575.42 
MHz. There is at least one persistent feature near 1581 
MHz, and perhaps more than one, but these are mixed 
among ripply interference that clearly fluctuates. It is 
therefore not possible to reliably estimate any of the 
OHM properties. 

IRAS F23050+0359: This OHM has no published 
spectrum , or propert i es, bu t was listed as an OH mega- 
maser bv lBaan et al.l ((1998). It is located at a redshift of 
0.0474, putting the OH lines at a frequency near a strong 
source of interference, likely ringing from GPS L2, that 
contaminates the low frequency side of the spectrum. 
Even so, there appears to be a detection of the 1667 MHz 
line, and possibly the 1665 MHz line. The Stokes V spec- 
trum has significant rippling, which prevents any possible 
Zeeman splitting detection or placement of useful limits 
on magnetic fields in this O HM. 

IRAS F23365+3604: IBottinelli etHI (fl990h re- 
ported this OHM to have an isotropic luminosity of 
log(LoH/-ko) = 2.44, but did not pub lish any othe r prop - 
erties of the source. Subsequently, iBaan et ail (1992) 
cited the upper limit on the isotropic luminosity of the 
1667 MHz line to be log(LoH/i©) < 1-93, assuming 
a 100 krns" 1 linewidth. iDarling & Giovanellil (|2002al ) 
considered the source suspect. The source is at a red- 
shift of 0.064531, which puts the 1667 MHz OH line at 
a frequency of 1566.3 MHz. The composite spectrum 
created by taking the mean of individual 1-second spec- 
tra shows significant broadband interference. Median 
combination of individual spectra to create a compos- 
ite removes most features, leaving what appears to be a 
reasonably flat spectrum with broad, low SNR emission 
from 1566.5-1568 MHz, plausibly close to the expected 
frequency of the OH line, and would correspond to an 
isotropic luminosity of log(ioH/£©) = 2.56. Spectra 
of the 1612 MHz line of IRAS F01417+1651, another 
source observed in this survey, spanned the frequency 
range 1566.1-1572.3 MHz. Examination of these spectra 
showed features of comparable shape and intensity to 
those that were appeared in the median composite spec- 
trum for IRAS F23365+3604, strongly suggesting that 
the apparent emission is actually interference. 

5.1.5. Stokes V non- detections 

For the sources described in this section, maser emis- 
sion was detected, but due to low signal-to-noise in 
the Stokes / and Stokes V spectra, it was not pos- 
sible to identify any features in the Stokes V as be- 
ing associated with features in the Stokes / spectrum. 
For the majority of sources, we do not present the 
spectrum, instead providing a reference where a pub- 
lished spectrum can be found. Exceptions are made for 
three sources, IRAS F08071+0509, IRAS F10339+1548, 
and IRAS F23028+0725. For IRAS F08071+0509, no 
Stokes / spectrum exists in the literature. In the case 
of IRAS F10339+1548, our data provided an improve- 
ment over the existing Stokes / data, allowing determi- 
nation of the hyperfine ratio, which is defined as the 
ratio of the flux in the 1667 MHz main line to the 
flux in the 1665 MHz main line, or Rh — F1667/F1665. 
IRAS F23028+0725 is included as an example of a source 
where our fitting routine reported marginally significant 
detection of Zeeman splitting that we do not find credi- 
ble. 



Where possible, we estimate an upper limit to the 
magnetic field str ength associated with the masing 
regions. Following iTroland fc Heilesl (|1982| ). the upper 
limit is taken to be the fitted field strength plus three 
times the error associated with the fitted field. Most 
OHMs are fit with multiple Gaussian components; for 
sources where an upper limit is reported, the value 
is that of the Gaussian component with the lowest 
associated error that appeared to be real. This is based 
on the shape and strength of features in the Stokes I 
spectrum, and the structure and amplitude of noise in 
the Stokes V spectrum. For example, limits associated 
with Gaussians fit to wide features in spectra are not 
reported, as they are likely bandpass features. We 
consider only one digit in these upper limits to be 
significant, but reported two digits for sources with 
upper limits between 10-20 mG. 

IRAS F01562+2528: This OHM, reported in DG02, 
has three broad, blended features, one of which is at 
the expected frequency of the redshifted 1665 MHz line. 
Despite the well-detected signal in Stokes /, there are no 
features in the Stokes V spectrum that look promising. 
Given the absence of any narrow structures, it is not 
possible to constrain the magnetic field in any meaningful 
way. 

IRAS 03521+0028: DG02 discovered this OHM, 
noting the excellent agreement between the peak of the 
1667 MHz emission and the published optical redshift, 
as well as marginal detection of the 1665 MHz line. Our 
spectrum has an additional narrow feature at 1448.3 
MHz not seen in the DG02 observations. Examination 
of multiple spectra for different objects that have obser- 
vations in this frequency range show that this new, nar- 
row feature is interference. No features in the Stokes V 
spectrum can be identified as associated with the astro- 
physical Stokes / emission, and the lack of bright, narrow 
emission prevents placing a useful limit on the magnetic 
field in this OHM. 

IRAS F03566+1647: This OHM was reported by 
DG02, and features broad, low amplitude emission. The 
Stokes V spectrum has no apparent features, and the 
broad Stokes I features do not provide useful limits on 
the magnetic fields in this source. 

IRAS F04121+0223: The spectrum for this OHM 
discovered by DG01 contains one narrow component cen- 
tered at 1485.20 MHz and one broad feature centered at 
1486.04 MHz. The Stokes V spectrum does not show 
any features associated with the Stokes / emission, and 
we place an upper limit of 40 mG on the magnitude of 
magnetic fields associated with the narrow emission in 
this OHM. 

IRAS 06487+2208: DG00 reported this OHM, fea- 
turing two main 1667 MHz peaks centered at 1457.7 MHz 
and 1458.3 MHz, and weak, corresponding 1665 MHz 
emission. There is a great deal of structure in the Stokes 
/ spectrum, requiring a total of nine Gaussians to pro- 
duce a good fit. The Stokes V spectrum features 1 mJy 
ripples, but no apparent Zeeman splitting signature. We 
place an upper limit of 10 mG on the magnetic field as- 
sociated with the narrow emission at 1457.7 MHz as well 
as the broader emission centered at 1458.3 MHz. 

IRAS 07163+0817: This OHM discovered by DG01 
features three distinct 1667 MHz peaks, one of which is 
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centered at 1501.3 MHz and is about twice as bright as 
the other two peaks. The Stokes V spectrum does not 
have any notable features, and we place an upper limit of 
20 mG on the magnitude of a magnetic field associated 
with the narrow line at 1501.3 MHz. 

IRAS 07572+0533: Only one broad feature can be 
distinguished for this OHM discovered by DG01. They 
noted an apparent second OH line offset blueward from 
the main feature by 400 km s~\ but we do not confi- 
dently detect this. Without any narrower lines, it is not 
possible to place useful limits on the magnetic fields in 
this source. 

IRA S F07 556+2859: This OHM was discovered by 
IWillett] (|2012() . and will be discussed in more detail in a 
future work by Willett, Darling, Kent, & Braatz. The 
Stokes I spectrum features broad, blended 1667 MHz 
emission and weak 1665 MHz emission. Without any 
distinct narrow emission, we do not place any limits on 
magnetic fields in this OHM. 

IRAS F08071+0509: IBottinelli et~aTl (fl989l) re- 
ported this OHM with log(L O H/£0) = 2.25, but did 
not p ublish a spe c trum or the flux of the 1667 MHz 
line. iBaan et al.l (|1998f ) classified it as a composite 
AGN/starburst in a survey to optically classify mega- 
maser galaxies. We clearly detect the maser, but there 
are no circu lar polarization features, as can be seen in 
Figure 5(a) The flux of the 1667 MHz line is 35 mJy, 
and the isotropic luminosity is log(LoH/ ^o) = 2.12. For 
their likely choice of cosmology, this would correspond 
to log(LoH/i©) = 2.03, which is roughly 70% less lumi- 
nous than they reported. It is not possible to identify 
Zeeman splitting in the Stokes V spectrum. There is a 
large dip in the Stokes V spectrum at the location of 
emission, but similar dips appear elsewhere in the spec- 
trum. The upper limit on the magnetic field associated 
with Gaussian 4, which is a bright, relatively narrow, 
feature located near the peak of emission, is 20 mG. The 
rest of the parameters used to fit emission are shown in 
Tabled 

IRAS F08201+2801: DG01 discovered this OHM, 
which features broad 1667 MHz emission, and low 
flux emission redward of the main 1667 MHz features, 
which DG01 suggested could be either 1665 MHz emis- 
sion or high velocity 1667 MHz emission. Overall, it 
has rich enough structure that ten Gaussian compo- 
nents are needed to produce a reasonable fit to the full 
1665/1667 MHz emission. The Stokes V spectrum has 
little in the way of structured noise, and shows no com- 
pelling features. We place an upper limit of 20 mG on 
fields in this OHM. 

IRAS F08279+0956: This OHM was reported by 
DG01. It features two broad, blended 1667 MHz features 
in the Stokes / spectrum, but at such low signal-to-noise 
that even with minimal structured noise in the Stokes V 
spectrum, it is not possible to place any useful constraint 
on magnetic fields associated with this OHM. 

IRAS F08449+2332: The maser emission of this 
OHM discovered by DG01 overlaps with the same 1448.3 
MHz interference that appeared in the spectrum of IRAS 
F03521+0028. As such, no magnetic field is detected, nor 
is it possible to place any useful limit on magnetic field 
properties. 

IRAS F08474+1813: DG01 discovered this OHM, 



which has a roughly 4 mJy peak flux density, and three 
broad emission features that are distinguishable. The 
low flux density and broad lines means the Stokes V ob- 
servations do not provide any useful limits on magnetic 
fields in this OHM. 

IRAS F09531+1430: This OHM was discovered by 
DG01, and the spectrum has two moderately narrow, 
well separated peaks atop a broad region of emission. 
Though no field is detected for either of the narrow fea- 
tures, we may place an upper limit of 30 mG on magnetic 
fields associated with either peak. 

IRAS F09539+0857: This OHM, discovered by 
DG01, has broad, blended 1665/1667 MHz emission and 
rich structure. We use nine Gaussian components to the 
total 1665/1667 MHz emission, including multiple nar- 
rower features, and achieve an adequate fit. The Stokes 
V spectrum has reasonably well behaved noise, but there 
are no suggestions of Zeeman splitting. The upper limit 
on the magnetic field associated with multiple narrow 
Gaussians fit to peaks in the emission is 20 mG. 

IRAS F10035+2740: There are two broad, low 
signal-to-noise features in the spectrum of this OHM, 
which was discovered by DG02. These broad features do 
not provide a useful limit on the magnetic fields in the 
masing region of this galaxy. 

IRAS F10339+1548: The OHM, discovered by 
DG01, is clearly detected at 1667 MHz. Weak 1665 MHz 
emis sion is also visible in our spectrum, shown in Figure 
5(b) [ while standing waves frustrated the detection of the 



1665 MHz emission in the DG01 observations. We fit one 
Gaussian to the 1665 MHz emission, and one broad and 
one narrow Gaussian to the 1667 MHz emission, which 
are shown in Table [3l From the Gaussian, we compute a 
hyperfine ratio Rh = 21. The peak of emission is fit by 
Gaussian 2, and with well behaved noise in the Stokes 
V spectrum, we can place an upper limit of 20 mG on 
magnetic fields in this source. 

IRAS F11028+3130: The spectrum for this OHM, 
reported by DG01, has two broad features, one corre- 
sponding to the 1667 MHz line and one to the 1665 MHz 
line, and a low amplitude 1667 MHz feature on the red- 
ward side of the broad 1667 MHz emission. Given the 
low snr, it is not possible to place any useful limit on the 
magnetic field strength for this source. 

IRAS F11180+1623: DG02 discovered this OHM, 
which has only one easily identified emission component 
with a peak flux density of roughly 4 mJy. With such a 
weak signal, no useful limit may be placed on the pres- 
ence of magnetic fields in this galaxy. 

IRAS F11524+1058: This OHM was discovered by 
DG01. It has broad emission, which we fit with four 
reasonably broad Gaussians. The lack of distinguishable 
narrow spectral features prevents placing useful limits on 
the magnetic field in this galaxy. 

IRAS F12005+0009: Though this OHM, discovered 
by DG02, has a peak flux density of only 9 mJy, it fea- 
tures multiple distinguishable emission components, in- 
cluding one relatively narrow peak centered on a broad 
hump. The fits to the Stokes V spectrum yields a 3c 
claimed detection for a 21 mG field associated with the 
narrow peak. Unfortunately, the Stokes V spectrum for 
this source includes many narrow ripples of the same 
characteristic size as the narrow feature, and shifting the 
Stokes / and V spectra relative to each other reveals 
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TABLE 2 

IRAS F08071+0509 Gaussian Fit Parameters 



S v Au vq Bu 

Gaussian (mjy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



2.74 ±0.20 1582.7638 ± 0.0143 0.4223 ± 0.0377 16023.2 -89.27 ±29.37 

1 3.03 ±0.21 1583.8481 ± 0.0130 0.3814 ± 0.0338 15807.0 -38.07 ± 25.31 

2 2.90 ±0.91 1584.5591 ± 0.0023 0.0902 ± 0.0082 15665.4 - 8.82 ± 13.27 

3 27.84 ±0.80 1584.6262 ± 0.0067 0.4756 ± 0.0243 15652.1 6.26 ± 3.23 

4 10.49 ±0.80 1584.7340 ± 0.0024 0.1168 ± 0.0082 15630.6 -9.77 ±4.19 

5 4.10 ±0.32 1585.2492 ± 0.0058 0.1532 ± 0.0142 15528.1 23.32 ± 11.83 



TABLE 3 

IRAS F10339±1548 Gaussian Fit Parameters 





S 


V 


Av 


Vq 


B \\ 


Gaussian 


(mjy) 


(MHz) 


(MHz) 


(km s — 1 ) 


(mG) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 





2.39 ± 0.14 


1391.4459 ±0.0198 


0.7723 ± 0.0550 


59446.6 


-136.05 ±49.28 


1 


6.69 ±0.41 


1393.1475 ±0.0101 


0.5373 ± 0.0203 


59007.8 


- 50.12 ± 15.57 


2 


11.26 ± 0.45 


1393.2540 ± 0.0025 


0.1841 ±0.0079 


58980.4 


- 2.35 ± 5.41 



many shifts that produce a claimed magnetic field de- 
tection. We place an upper limit of 40 mG on magnetic 
fields associated with masing clouds in this galaxy. 

IRAS F12018+1941: Only one broad feature may 
be distinguished i n the sp ectrum of this OHM discovered 
by iMartin et al.l (|1988bf) . It is not possible to place a 
useful limit on the magnetic field in this OHM. 

IRAS F12162+1047: DG02 reported this OHM, 
which has broad, low flux density emission. The Stokes 
/ spectrum may be adequately fit with two broad Gaus- 
sians, but does not allow any useful limit to be placed on 
magnetic fields in this gal axy. 

IRAS F12243-0036: IMartin et ail (|1988bD discov- 
ered this OH kilomaser, which has apparent absorption 
bracketing the maser emission. The overall profile is com- 
plicated, and may only be adequately fit. Nevertheless, 
there are a few narrow emission features in the spectrum, 
and one that is aligned with a feature in the Stokes V 
spectrum. The fit to Stokes V suggests a field detection, 
but shifts between Stokes / and V also produce claimed 
fields. We place an upper limit of 20 mG on the magnetic 
field associated with this maser. 

IRAS F12549+2403: DG02 discovered this OHM, 
which has two broad, 1667 MHz emission components 
and marginally detected 1665 MHz emission in the Stokes 
/ spectrum. No useful limit may be placed on the mag- 
netic field in this galaxy. 

IRAS F13126+2453: This s ource is classified as 
an OH absorber and kilomaser. iDarling fc Giovanellil 
(|2002af) reported that no published spectrum exis t s, but 
a published spectrum is available in lSchmelz et al.l (1986) 
under the designation IC 860, so the spectrum is not pro- 
vided again here. The absorption in both the 1665 and 
1667 MHz lines is clearly detected, as well as emission at 
the wings of the absorption features. For the absorption 
alone, the hyperfine ratio is Rr = 1-27 ± 0.18, while for 
the emission, we find a hyperfine ratio of Rh = 1.3±0.36, 
where the errors are based on reported errors to Gaus- 
sian fits to the spectra. Both are consistent with OH in 
local thermodynamic equilibrium. 



The Stokes V spectrum features bandpass structure 
on scales comparable in width to features in the Stokes 
/ spectrum, leading to claimed fields with magnitudes of 
~30-40 mG associated with the absorption features. We 
do not consider these claims to be reliable. The noise in 
the vicinity of the emission components has less structure 
than in the regions of absorption, but there is no sign 
of Zeeman splitting. In the emission regions, we place 
an upper limit of 30 mG on the magnitude of magnetic 
fields. 

IRAS F13218+0552: Among broad emission in the 
spectrum of this OHM discovered by DG02, there is one 
narrow feature at 1384.8 MHz. The Stokes V spectrum 
has well-behaved noise in the vicinity of the narrow emis- 
sion seen in Stokes J, but there is no sign of Zeeman split- 
ting. We place an upper limit of 50 mG on the magnetic 
field associated with the masing region in this source. 

IRAS F14043+0624: This OHM, discovered by 
DG02, has an asymmetric 1667 MHz emission profile best 
fit by one narrow and one moderately wide Gaussian, 
along with one broad, weak 1665 MHz line. No magnetic 
field is detected, and we place an upper limit of 40 mG 
on magnetic fields in this galaxy. 

IRAS F14059+2000: DG02 discovered this OHM, 
which has bright 1667 MHz emission that is well fit 
by one very broad Gaussian and one moderately broad 
Gaussian, and 1665 MHz emission that blends with the 
red wing of 1667 MHz emission. No Zeeman splitting is 
apparent in the Stokes V spectrum, and we place an up- 
per limit of 20 mG on the magnitude of magnetic fields 
associated with this OHM. 

IRAS F14553+1245: There is one relatively narrow, 
low SNR feature in the spectrum of this OHM. It was 
discovered by DG02, and like them, we do not see a clear 
1665 MHz component in the spectrum. The Stokes V 
spectrum is nearly featureless, and we place an upper 
limit of 60 mG on magnetic fields associated with this 
OHM. 

IRAS F14586+1432: This OHM, discovered by 
DG02, has a rich, blended spectrum, spread over more 
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than 5 MHz. Despite this richness, there are few identi- 
fiable narrow features in the spectrum. Fits to Stokes V 
provide claimed magnetic fields, but so too do fits when 
the Stokes V is shifted relative to Stokes /, thanks to a 
number of ripples in the Stokes V spectrum. This is re- 
garded as a non-detection, and it is difficult to provide a 
meaningful upper limit on magnetic fields in this ga laxy. 

IRAS F15107+0724: Bottinelli et al.l (fl986l) dis- 
covered th i s OH M, a nd a spectrum is available in 
iBaan et al.l (|1987l ) and iMartin et al.1 (|1988af ). The 1667 
MHz emission has a dual-peaked structure, which is 
closely matched by the 1665 MHz emission. A total of 
six Gaussian components are needed to fit the Stokes / 
spectrum, four of which are narrow. While no Zeeman 
splitting is seen in the Stokes V spectrum, each of the 
narrow components is consistent with an upper limit of 
12-15 mG to the magnitude of magnetic fields in this 
OHM. 

IRAS F16100+2527: DG01 discovered this OHM, 
which has emission made up of a broad 1667 MHz line 
with a marginally distinct peak, and a similarly broad 
and well separated 1665 MHz line. The 1667 MHz emis- 
sion can be reasonably fit with one wide and one narrow 
Gaussian. The nearly 6 mJy narrow feature provides a 
40 mG upper limit on the magnitude of magnetic fields 
in this OHM. 

IRAS F16300+1558: This OHM was discovered by 
DG00, but RFI plagued their observations. Our spectra 
combined by taking the median do not have serious nar- 
row RFI, but do feature significant bandpass structure. 
Like DG00, we do not clearly detect the 1665 MHz line. 
Structure of unclear origin in both the Stokes / and V 
spectra prevents placing any useful limits on magnetic 
fields in this OHM. 

IRAS F17161+2006: The 1667 MHz emission from 
this source, discovered by DG02, has a peak flux density 
of less than 10 mJy. Four Gaussians are required to rea- 
sonably fit the profile, three for the 1667 MHz emission 
and one for the broad 1665 MHz emission. One of the 
1667 MHz components is reasonably narrow, and pro- 
vides an upper limit of 60 mG on the magnitude of the 
magnetic field associate d with this OHM . 

IRAS F17207-0014: iBottinelli et~aT1 (|1985D reported 
this OHM, which has a spectrum with strong interfer- 
ence from GLONAS. Despite the interference, the maser 
is still clearly detected, as interference is negligible at fre- 
quencies lower than about 1601 MHz, based on inspec- 
tion of individual spectra used to create the composite 
spectrum. Comparison of previousl y published spectra 
for the OHM in IRA S F17207-0014 ([Martin et al.lH98l 
iMomiian et al.l I2006D , as well as with spectra of other 
sourcing overlapping in frequency, confirms that the fea- 
tures below 1600 MHz are maser emission rather than 
interference. The Stokes / profile is bright, but there 
are no distinct peaks in the emission. Moreover, the 
Stokes V spectrum has significant broadband structure. 
While interference below 1601 MHz is mild compared to 
that above 1601 MHz, it still makes it difficult to detect 
Zeeman splitting in this source or place useful limits. 

IRAS F17539+2935: This OHM, discovered by 
DG00, is only marginally detected, with a peak flux den- 
sity of just over 1 mJy. The SNR is insufficient to detect 
magnetic fields or place any useful upper limits. 

IRAS 20248+1734: While this OHM, discovered by 



DG00, is clearly detected, the spectrum for this OHM 
features standing waves that complicate fitting and iden- 
tification of real emission, and it is not possible to detect 
Zeeman splitting or place useful limits on magnetic fields 
in this galaxy. 

IRAS 20286+1846: DG00 discovered this OHM, 
which has extremely broad, blended 1665/1667 MHz 
masing lines. To fit the total profile adequately requires 
nine components, three of which are relatively narrow. 
No Zeeman splitting is detected for any of the three nar- 
row components, but we place an upper limit of 20 mG 
on magnetic fields in this OHM. 

IRAS 20450+2140: This OHM was reported in 
DG00. The Stokes / spectrum shows two blended 
1667 MHz emission components, adequately fit by two 
Gaussians. The width of the fitted Gaussians is too 
large to place any useful limits on magnetic fields in this 
galaxy. 

IRAS 21077+3358: DG00 reported this OHM, 
which has broad 1667 MHz that possibly blends with its 
1665 MHz component. It is not possible to clearly dis- 
tinguish any narrow features in the spectrum, or place 
any useful limits on magnetic fields in this galaxy. 

IRAS 21272+2514: The spectrum of this OHM, dis- 
covered by DGOO, is rich in structure, and very broad. It 
is not possible to clearly distinguish the 1665 MHz and 
1667 MHz lines. The full spectrum is reasonably well-fit 
by ten Gaussians, six of which are moderately narrow. 
No magnetic fields are detected, and we place an upper 
limit of 20 mG on magnetic fields in the masing region 
in this galaxy. 

IRAS 22055+3024: The spectrum of this OHM, re- 
ported in DG01, has moderately broad 1667 MHz emis- 
sion and a distinct 1665 MHz line. The 1667 MHz fea- 
tures are well fit with a wide base and two narrow Gaus- 
sians. There is some rippling structure in the Stokes V 
in the vicinity of the narrow lines, but the noise is well 
enough behaved to provide a 30 mG upper limit on mag- 
netic fields in this OHM. 

IRAS F22116+0437: This OHM was detected by 
DGOO. Though it is re-detected here, it is insufficiently 
bright to place any useful limit on magnetic fields. 

IRAS F23019+3405: Emission from this OHM, re- 
ported in DG01, is quite narrow, with a FWHM of 
roughly 60 kms 1 . It nevertheless appears to have a 
bi-peaked structure, though higher SNR spectra would 
be required to confirm this. No Zeeman splitting is ap- 
parent in the Stokes V spectrum, but we are able to place 
an upper limit of 30 mG on magnetic fields in this OHM. 

IRAS F23028+0725: DG01 reported this OHM, 
which has a broad 1667 MHz emission that blends with a 
clear 1665 MHz component. Five Gaussians are needed 
to adequately fit t he co mbined 1665/1667 MHz emission, 



shown in Figure 5(c) Gaussian components 1 and 2 
are moderately narrow, and produce fits with reported 
~ 3<r magnetic fields having magnitudes of ~ 40 mG. 
The spectrum is not convincing to the eye, however, and 
shifting the Stokes / and V relative to each other con- 
firms this, as non-zero shifts produce stronger fits to the 
Stokes V spectrum, as seen in Figure [TJ For this reason, 
it is difficult to place useful limits on magnetic fields in 
this galaxy. 

IRAS F23129+2548: Emission in this OHM, re- 
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ported in DG01, is broad, and the 1665 MHz and 
1667 MHz lines are blended. The total spectrum is well 
fit with six Gaussians, though only one is narrow. Zee- 
man splitting is not observed, but noise in the Stokes 

V spectrum is well enough behaved to provide an upper 
limit of 50 mG on magnetic fields in this OHM. 

IRAS F23199+0123: This OHM was reported by 
DG01. The 1667 MHz emission is broad, and has a peak 
flux density of only 3 mJy It is not clear if 1665 MHz 
emission is detected. Given the low flux density, no 
meaningful limit may be placed on magnetic fields in 
this OHM. 

IRAS F23234+0946: DG01 discovered this OHM. It 
has broad 1667 MHz emission with an apparent blue tail 
and 1665 MHz line. Four total Gaussians adequately fit 
the spectrum, one of which is narrow and near the peak 
of 1667 MHz emission. No Zeeman splitting is detected, 
but the narrow feature provides an upper limit of 50 mG 
on the magnitude of magnetic fields in this OHM. 

5.1.6. Marginal Stokes V detections 

The three sources presented in this section, 
IRAS F15224+1033, IRAS F15587+1609, and 
IRAS F20550+1655, all have features in their Stokes 

V spectra consistent with Zeeman splitting. They also 
have Stokes V structure that does not appear to be 
associated with any maser emission. For this reason, 
they are considered only "marginal" detections. In 
Table HJ Table and Table d the fitted features that 
are consistent with detections of Zeeman splitting are 
highlighted with bold text. 

IRAS F15224+1033: This OHM, discovered by 
DG01, has a rich, wide Stokes / spectrum, and a fairly 
noisy Stokes V spectrum, shown in Figure 5(d) Any 
1665 MHz emission that might be present is blended with 
the broad 1667 MHz emission. The features of the OHM 
have not changed significantly since DG01. The complex- 
ity requires nine Gaussians to adequately fit, and even 
then, the residuals show structure that coincides with the 
two brightest features, which are fit by Gaussian compo- 
nents 3, 4, and 6. The full fits are shown in Table |4j 
At the same frequency as Gaussian 6, there is a weak 
feature in the Stokes V spectrum, with a peak-to-trough 
amplitude only marginally larger than that of others in 
the spectrum. The fits regard this as a 3.9<r detection of 
a -8.3 mG magnetic field. There is also a marginally sig- 
nificant fit associated with Gaussian 4, with a reported 
field of -26.7±8.5 mG, but the fitted structure in Stokes V 
does not appear markedly different to eye than other fea- 
tures in the spectrum. Altogether, the structure in the 
Stokes V spectrum that is apparently unassociated with 
features in the Stokes I spectrum, and the marginal ap- 
pearance of the fitted features, limits our confidence in 
the two reported fields. 

IRAS F15587+1609: The spectrum of this OHM, 
reported in DG01, has a great deal of structure, with 
thr ee dis tinct peaks, and many narrow features. Fig- 
ure 5(e) shows our fits, and Table [5] lists the fits, which 
use eleven total Gaussians for the 1667 MHz emission, 
and two more for the 1665 MHz emission. In fitting 
the Stokes V spectrum, there are three reported mag- 
netic fields, associated with Gaussians 6, 7, and 8, and 
each is just above 3cr claimed detections. These are 
three of the four brightest Gaussian components used, 



have similar reported values of 6.1 ± 2.0, 7.5 ± 1.9, 
and 12.2 ± 3.5, respectively, and are contained within 
a roughly 100 km s _1 range of velocities. 

We regard this claimed magnetic field detection as 
marginal primarily because the fitted features do not look 
particularly compelling to the eye. Though there are no- 
ticeable dips in the Stokes V spectrum near the same 
frequency as the brightest emission seen in the Stokes V, 
there are no features strongly resembling Zeeman split- 
ting. Given the complex structure of the spectrum, shift- 
ing the Stokes / and V spectra relative to each other does 
not provide significant insight. The fits are of highest 
quality for no shift, but are of only slightly lower quality 
at a shift of 0.12 MHz, which is roughly the average of 
the frequency separation between Gaussian 5, 6, 7, and 8. 
Shifting the Stokes I and V by roughly 0.5 MHz relative 
to each other, corresponding to the separation between 
the central peak and the redward peak, also produces 
claimed fields, though they are also of lower quality than 
when no shift is applied. Taken together, we are not 
confident in these reported fields. 

IRAS F20550+ 1655: This sou r ce wa s reported in 
an IAU Circular by IBottinelli et al.l (119861 ) , and a pub- 
lished spectrum is availab le inlBaan et al.l (119891) . under 
the designation II Zw 96. IBottinelli et al.l (|1986l ) cite an 
isotropic luminosity of log(LoH/ Lq) — 2.27, and later 
sources use log(LoH/ '£©) = 2.11, in agreement with the 
value of log(LoH/£©) = 2. 13 that we find. The spec- 



trum, shown in Figure 5(f) suffers from interference in 
the 1602-1608.5 MHz part of the spectrum that is caused 
by the GLONAS satellite, including apparent absorption 
that just overlaps with redward side of the 1667 MHz line, 
rcdshifted to 1609 MHz. The interference is strongly lin- 
early polarized, and also produces features in the Stokes 
V spectrum. It is not clear if this interference appears 
beyond 1608.5 MHz in either the Stokes I or V spectra. 
Using eleven Gaussian components, provided in Table [6j 
to fit the emission blueward of 1608.5 MHz does an ad- 
equate job, and produces four magnetic field fits greater 
than 3<T. Two of these are to wide components, and 
are not considered believable, but Gaussians 2 and 5 are 
narrow, and produce claimed fields of —18.5 ± 2.4 mG 
and — 17.4 ± 2.9 mG, respectively, and also appear plau- 
sible examining the fits visually. The two Stokes V fea- 
tures that are fit are the largest amplitude, from peak- 
to-trough, of any blueward of 1608.5 MHz, but there are 
two other dips that are not well accounted for. The inter- 
ference at 1608 MHz is also shown, and is of comparable 
amplitude to the fitted features. Though it has a differ- 
ent shape, it casts doubt on the validity of the detection. 
Examining the variance of each channel of the Stokes V 
spectrum shows a local maximum at 1608 MHz, corre- 
sponding to the interference, but has no apparent struc- 
ture blueward of 1608.5 MHz, which suggests the fitted 
features in the Stokes V spectrum are real. Given the 
conflicting evidence, we consider this a plausible mag- 
netic field detection. 

5.1.7. Stokes V detections 

The following sources all exhibit features in their 
Stokes V spectra that are consistent with Zeeman split- 
ting of OH maser lines. In the tables presented for 
sources with detections, the features that are considered 
detections of Zeeman splitting are highlighted with bold 
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TABLE 4 

IRAS F15224+1033 Gaussian Fit Parameters 







S 






V 






At/ 


Vq 




By 


Gaussian 




(mjy) 






(MHz) 






(MHz) 


(km s~ 




(mG) 


I 1 ) 




( 2 ) 






(3J 






(4) 


( 5 ) 




W 





1 


35 ±0 


18 


1467 


4444 ± 


0212 





3575 ± 0.0593 


40841 


7 


-179.44 ±43.83 


1 


1 


11±0 


23 


1468 


0327 ± 


0197 





2063 ± 0.0520 


40705 


2 


52.26 ± 40.79 


2 


3 


25 ±0 


19 


1468 


7777 ± 


0272 





4532 ± 0.0595 


40532 


5 


- 40.08 ± 23.32 


3 


5 


36 ±0 


69 


1469 


1573 ± 


0293 





2430 ± 0.0619 


40444 


5 


- 14.76 ± 10.59 


4 


4 


75 ±1 


78 


1469 


3149 ± 


0118 





1446 ± 0.0255 


40408 





-26.62 ±8.45 


5 


8 


30 ±0 


32 


1469 


7945 ± 


0010 





8947 ± 0.0029 


40297 





-0.09 ±12.39 


6 


15 


76 ±0 


22 


1469 


7968 ± 


0094 





1116 ±0.0338 


40296 


5 


-8.27±2.13 


7 


4 


34 ±0 


29 


1470 


5376 ± 


0686 


6 


3929 ± 0.4799 


40125 


2 


86.82 ± 66.19 


8 


6 


07 ±0 


20 


1470 


8441 ± 


0059 





3314 ±0.0156 


40054 


3 


14.55 ± 9.88 


9 


2 


09 ±0 


19 


1471 


3402 ± 


0175 





3661 ± 0.0478 


39939 


7 


- 8.39 ± 30.00 



TABLE 5 

IRAS F15587±1609 Gaussian Fit Parameters 







S 






V 






At/ 


VQ 






Gaussian 




(mjy) 






(MHz) 






(MHz) 


(km s~ 


') 


(mG) 


(1) 




(2) 






(3) 






(4) 


(5) 




(6) 





1 


91 ±0 


28 


1464 


5345 ± 


0130 





2546 ± 0.0426 


41117 


9 


24.22 ± 16.49 


1 


6 


23 ±0 


14 


1465 


0808 ± 


0129 





8858 ± 0.0418 


40990 


7 


-14.98 ±9.46 


2 


7 


08 ±0 


53 


1466 


3288 ± 


0038 





0686 ± 0.0076 


41100 


8 


0.22 ±4.21 


3 


7 


75 ±0 


44 


1466 


4171 ± 


0039 





0862 ± 0.0100 


41080 


3 


0.84 ±4.32 


4 


15 


30 ±0 


61 


1466 


5661 ± 


0337 





7174 ± 0.0456 


41045 


7 


- 2.04 ± 6.34 


5 


7 


06 ±0 


87 


1466 


6249 ± 


0045 





0843 ±0.0122 


41032 





12.01 ±4.90 


6 


18 


78 ±1 


10 


1466 


7254 ± 


0017 





0847 ± 0.0053 


41008 


7 


6.12 ±2.01 


7 


25 


75 ±3 


70 


1466 


8546 ± 


0036 





1567 ±0.0125 


40978 


6 


7.51 ± 1.94 


8 


20 


80 ±1 


71 


1467 


0327 ± 


0110 





2574 ± 0.0458 


40937 


3 


12.23 ±3.48 


9 


4 


19 ±0 


64 


1467 


2375 ± 


0026 





0482 ± 0.0085 


40889 


7 


-6.87 ±5.47 


10 


25 


91 ± 1 


29 


1467 


3137 ± 


0076 





2382 ± 0.0130 


40872 





0.60 ±2.57 


11 


6 


45 ±0 


17 


1467 


6886 ± 


0156 





5179 ± 0.0375 


40785 





16.38 ±13.72 


12 


2 


97 ±0 


13 


1468 


7666 ± 


0258 


1 


2354 ± 0.0881 


40535 





-65.71 ±40.91 



TABLE 6 

IRAS F20550±1655 Gaussian Fit Parameters 







S 






V 






Ai/ 






By 


Gaussian 




(mjy) 






(MHz) 






(MHz) 


(km s~ 


J ) 


(mG) 


(1) 




(2) 






(3) 






(4) 


(5) 




(6) 





7 


17 ± 1 


14 


1608 


6408 ± 


0014 





0180 ± 0.0034 


10942 


9 


-2.27 ±3.20 


1 


31 


72 ± 1 


10 


1608 


9151 ± 


0050 





2567 ± 0.0085 


10890 





-10.01 ±2.91 


2 


17 


19 ± 1 


06 


1608 


9286 ± 


0019 





0481 ± 0.0040 


10887 


1 


-18.53 ±2.39 


3 


15 


47 ± 1 


50 


1608 


9974 ± 


0024 





0693 ± 0.0095 


10874 


1 


- 2.73 ± 3.13 


4 


28 


65 ± 1 


32 


1609 


0765 ± 


0016 


1 


1363 ± 0.0044 


10858 


8 


-31.89 ± 7.26 


5 


10 


13 ±0 


91 


1609 


1010 ±0 


0049 





0276 ± 0.0139 


10854 


1 


-17.42 ±2.87 


6 


5 


91 ±0 


66 


1609 


1543 ± 


0064 





0635 ±0.0216 


10843 


8 


-15.02 ± 7.42 


7 


9 


12 ±0 


83 


1609 


2766 ± 


0018 





0493 ± 0.0056 


10820 


2 


-5.59 ±4.02 


8 


12 


72 ±0 


63 


1609 


4236 ± 


0017 





0810 ± 0.0054 


10791 


8 


-7.27 ±3.69 


9 


3 


20 ±0 


56 


1609 


6839 ± 


0067 





0820 ±0.0180 


10741 


6 


48.31 ± 14.76 


10 


3 


70 ±0 


68 


1609 


8747 ± 


0046 





0515 ±0.0115 


10704 


8 


-4.71 ± 10.14 



text. In some cases, however, it was not possible to re- 
liably estimate the magnetic field strengths associated 
with maser components. Of the new detections, none of 
them have any existing VLBI observations, so it is not 
possible to provide detailed comments on the structure 
of the magnetic fields in the OHMs. 

IRAS F02524+2046: DG02 discovered this very lu- 
minous OHM, which features three regions of narrow 
1667 MHz emission, and two distinct areas of 1665 MHz 
emission that corre spond to two of the 1667 MHz emis- 
sion regions. Figure 5(g) shows a Stokes V spectrum that 
is similarly rich, relative to the other Stokes V spectra 



of sources in this survey. There are three large ampli- 
tude features right next to one another, at frequencies 
corresponding to the brightest emission in the Stokes / 
spectrum. 

Fitting all of this structure was challenging. Given 
the relatively lower flux of the 1665 MHz emission, five 
Gaussian components were sufficient to produce resid- 
uals absent of any structure. Most structure in the 
1667 MHz emission in the Stokes I spectrum could be 
reasonably fit with nine or ten Gaussian components, 
but produced extremely poor fits to the most obvious 
features in Stokes V, and left residuals around some of 
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the narrower peaks. Ultimately, we used thirteen compo- 
nents for fitting 1667 MHz features, which are all listed 
in Table [3 and produced a fit to Stokes I that still has 
difficulty with the brightest feature in the spectrum, and 
only partially captures the complex Stokes V structure. 
In particular, a strong dip near 1412.1 MHz is completely 
missed by the fits, and no combination of Gaussians we 
attempted was able to simultaneously fit that feature and 
the other strong features in the Stokes V spectrum of this 
source. 

Nevertheless, the fields and errors associated with 
these fits yield multiple confident detections, as well as 
two fits that exceed 3er but in which we have no confi- 
dence. The claimed, but ultimately dubious, fields are 
those associated with Gaussians 2 and 10. Each of these 
two Gaussians is fitting broad emission at a velocity of 
roughly 54,100-55,250 km s~\ Gaussian 2 for 1665 MHz 
emission and Gaussian 10 for 1667 MHz emission. We 
generally place no confidence in such wide features, as 
discussed in Section l4.2.1l 

Aside from these reported fields, there are four strong 
fields associated with Gaussians 8, 9, 12, and 13. All 
are narrow Gaussians fit to the bright, central 1667 MHz 
emission, and each gives a strong positive field, rang- 
ing from 12.3±1.8 mG for Gaussian 9 to 23.9±5.8 mG 
for Gaussian 13. While this result is dependent upon 
the input parameters of the fit, different combinations of 
Gaussians consistently reported strong, positive fields as- 
sociated with the masing regions that produced the cen- 
tral 1667 MHz emission. A weak 1665 MHz feature, fit 
by Gaussian 3, provides further confirmation of a strong, 
positive field. On its own, this weak component and 
< 3<7 fitted field would not be considered significant, but 
it has a similar velocity and width as some of the signif- 
icant 1667 MHz features. The component is best fit by 
a field of 15.5 ± 5.2 mG, right in the same range of mag- 
netic field strengths measured by the 1667 MHz lines. 
VLBI observations of this OHM do not, unfortunately, 
exist, which limits further interpretation of the structure 
in the masing region. 

IRAS F04 332+0209: This OH kilomaser was first 
reported by iMartin et al.l (|1989h with log(L H/£©) = 
0.44, but no spectrum or other properties of the maser 
were publ ished . We clearly detect this source, shown 
in Figure 5(h) It has a 1667 MHz peak flux density 



of 16 mjy, and an isotropic luminosity log(LoH/£©) 
0.29, assuming a distance of 51 Mpc. This i s slig htly 
lower than the value reported by [Martin et al.l (|1989l ). If 
they assumed a distance of 48 Mpc, our measured flux 
suggests possibility that the luminosity of the source has 
changed of order 40-50%. Given the unknown uncer- 
tainty in their luminosity estimate, uncertainty in their 
assumed distance, and roughly 10% uncertainty in ours, 
this evidence for variability is merely suggestive. 

As a kilomaser host, this galaxy is unlike the majority 
of the galaxies observed in this survey. It has an in- 
frared luminosity of 2.9 x10 10 Lq, well short of meeting 
the definition of a LIRG, which requi res log(LFm/£o) > 
11. Nevertheless, iBaan et al. (1998) characterized this 
galaxy as having a starburst nucleus in their optical clas- 
sification of maser hosts. 

The maser is also unusual in that the Stokes V spec- 
trum for this maser is unlike any other in this survey. 



There are three distinct features in the spectrum at the 
same frequency as emission in the Stokes / spectrum, and 
the features in Stokes V have the same width and magni- 
tude of the amplitude as those in the Stokes / spectrum. 
This suggests that the narrow splitting assumption is in- 
valid for this source; instead, the Zeeman splitting ex- 
ceeds the lincwidth, and the LHCP and RHCP compo- 
nents are distinct. In the absence of VLBI observations 
for this source, this makes deriving estimates of magnetic 
fields difficult. For this reason, we do not present a table 
showing fits, as we do for other sources. 

Nevertheless, we offer a plausible scenario, while ac- 
knowledging that this single dish spectrum provides in- 
complete evidence. The two reddest peaks in the Stokes 
/ spectrum, centered at 1647.43 MHz and 1647.53 MHz 
(labelled on the axis of the spectrum with a 0), have com- 
parable widths, and amplitudes that differ by a factor of 
roughly two. The 1647.43 MHz peak is the weaker of the 
two, and almost entirely RHCP, while the 1647.53 MHz 
peak is LHCP. If the two are paired components, then 
the difference between centers of the Gaussian compo- 
nents, Ai/, gives the total magnetic field 



D 



Aiv 
IT' 



C-i) 



where b is the same splitting coefficient described be- 
fore. The roughly 100 kHz separation would correspond 
to a 46.8 mG field. The Stokes V spectrum contains 
another large dip at the same frequency as the bluest 
feature in the Stokes / spectrum. This component is 
nearly all in the LHCP spectrum. If the component is 
associated with a field that is also strong, and positive, 
the corresponding RHCP emission would fall near the 
central, brightest emission from this maser. The Stokes 
V spectrum is relatively flat in the region of the brightest 
emission. This could result from overlap of the RHCP 
emission associated with the reddest feature and LHCP 
emission associated with the peak, along with asymmet- 
ric amplification of the LHCP and RHCP components 
in the source. Even if this interpretation is incorrect, 
we regard this as a clear detection of Stokes V features. 
Follow-up VLBI observations would be very interesting 
in better understanding this source. 

IRAS F09039+0503: DG01 reported this OHM, 
identifying weak 1665 MHz emission as well as a broad 
blueshifted 1667 MHz component at 37,300 km s _1 
in addition to the main features between 37,500- 
37,900 km s _1 . The 1665 MHz feature is apparentl y ab- 



sent in this more recent spectrum, shown in Figure 5(i) 
and the blueshifted 1667 MHz component is marginally 
detected. The central part of the spectrum features mul- 
tiple narrow peaks, ideal for Zeeman detections. The 
Stokes V spectrum shows two clear features, each asso- 
ciated with the bright, narrow peaks, and fit by Gaus- 
sians and 2. The fitted magnetic fields are each con- 
sidered significant by the fits, and none of the checks we 
tried suggested that the features were not real. The field 
reported for Gaussian is -16.1±2.7 mG, while Gaus- 
sian 2 has an associated field of -27.6±5.5 mG. The 
Stokes V feature fit by Gaussian 2 is clearly asymmet- 
ric, possibly as a result of the relatively wide splitting 
relative to the line width. None of the other fitted com- 
ponents, shown in Tabled are considered significant. Of 
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TABLE 7 

IRAS F02524+2046 Gaussian Fit Parameters 





S 


V 


At/ 


Vq 


By 


Gaussian 


(mjy) 


(MHz) 


(MHz) 


(km s _1 ) 


(mG) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


a 


11.37 ± 0.68 


1410.0094 ± 0.0014 


0.0910 ± 0.0053 


54300.9 


1.36 ± 2.41 


1 


6.87 ± 0.67 


1410.0398 ± 0.0078 


0.2969 ± 0.0259 


54293.2 


— 13.10 ± 7.58 


2 


11.07 ± 0.36 


1410.4608 ± 0.0044 


0.3152 ± 0.0122 


54187.5 


— 17.83 ± 4.49 


3 


3.63 ± 0.47 


1410.6113 ± 0.0030 


0.0510 ± 0.0083 


54149.8 


15.50 ± 5.19 


4 


7.40 ± 0.28 


1410.8600 ± 0.0276 


1.3718 ± 0.0806 


54087.4 


21.66 ± 13.77 


5 


17.43 ± 0.74 


1411.6473 ± 0.0056 


0.1795 ± 0.0080 


54305.7 


7.50 ± 3.93 


6 


7.10 ±4.52 


1411.8708 ± 0.0438 


0.2179 ± 0.1105 


54249.6 


- 1.15 ± 12.19 


7 


9.53 ±3.46 


1412.0160 ±0.0029 


0.0562 ± 0.0098 


54213.2 


- 6.73 ± 3.86 


8 


38.02 ±7.69 


1412.0955 ±0.0057 


0.1146 ±0.0185 


54193.3 


16.22 ± 1.61 


9 


32.04 ±8.95 


1412.2010 ±0.0089 


0.1133 ±0.0155 


54166.8 


12.27 ± 1.83 


10 


35.20 ±5.23 


1412.2861 ±0.0513 


0.4675 ± 0.0484 


54145.5 


-12.07 ±3.18 


11 


11.97 ±1.02 


1412.3222 ±0.0007 


0.0230 ± 0.0022 


54136.4 


-0.78 ±1.85 


12 


41.25 ±2.55 


1412.3425 ±0.0012 


0.0813 ±0.0027 


54131.4 


13.65 ± 1.07 


13 


5.56 ±0.52 


1412.5320 ±0.0021 


0.0538 ± 0.0061 


54083.9 


23.88 ±5.82 


14 


4.28 ±0.50 


1412.7663 ±0.0023 


0.0411 ±0.0060 


54025.2 


10.94 ± 6.63 


15 


12.57 ±0.48 


1412.9512 ±0.0012 


0.0832 ± 0.0038 


53978.9 


1.75 ± 3.33 


16 


11.75 ±0.38 


1412.9866 ±0.0055 


0.3496 ±0.0161 


53970.0 


-9.07 ±7.83 


17 


3.60 ±0.22 


1413.5314 ±0.0269 


0.6639 ± 0.0795 


53833.7 


75.98 ± 33.51 



a Velocities and magnetic fields for Gaussians 0-4 correspond to values for the 1665 MHz line. 



all sources in the survey with magnetic fields detections 
in which we are confident, the fields reported here for 
IRAS F09039+0503 are the strongest. 

IRAS F10378+1 108: This OHM was discovered by 
iKazes fc Baanl (fl99ll ). and re-observed by DG02. The 
spectrum, shown in Figure 5(j)| features extremely broad 
emission, with two distinct narrow features overlaid. The 
broad spectrum requires 12 Gaussians for a reasonable 
fit, listed in Table [9j though the residuals show that the 
fits near the two narrow features could still be improved. 
The Stokes V spectrum contains two suggestive features. 
One of these is well aligned with the peak in Stokes / fit 
by Gaussian 3, and the field associated with this feature 
is -13.9 ± 3.3 mG. The other apparent feature is just 
redward of the peak fit by Gaussian 7, and is fit in part by 
Gaussian 6, and a fitted field of -28.0 ± 7.0 mG associated 
with it. Given the width of Gaussian 6, we do not regard 
this reported field as being meaningful. 

IRAS F16255+2801: This OHM was discovered and 
characterized by DG01. Like DG01, we only detect 
1667 MHz emission. The Stokes / spectrum, shown in 



regard this as a clear detection of Zeeman splitting, we 
are less confident in the derived magnetic field strengths. 

IRAS 18368+3549: The emission from this OHM 
discovered by DG01 is extremely broad, causing any 
1665 MHz to be blended with the 1667 MHz emission, 



Figure 5(k) (with fits listed in Table [TUt . features four 
main overlapping peaks, and there is significant structure 
in the Stokes V spectrum at the same frequency. The 
Stokes / could be fit well with five Gaussian components, 
but the addition of one more low amplitude, narrow com- 
ponent on the red side of the spectrum significantly im- 
proved the fit to the Stokes V spectrum, and even then 
the Stokes V fit is still poor. There is clear, strong struc- 
ture, and the fields associated with Gaussian 1, 2, and 
3 are all significant: 7.0 ± 1.9 mG, 9.6±2.0 mG, and 
26.8±4.0 mG. Given the large field claimed for Gaussian 
3, and the relative narrowness of the lines, the narrow 
line assumption we use in fitting fields is perhaps pushed 
past its limit, as was the case for IRAS F04332+0209. We 
experimented with separately fitting the RCP and LCP, 
using only five Gaussian components, and deriving mag- 
netic field estimates using the frequency separation. The 
results are actually fairly similar, though the reported 
fields for each component are 20-30% smaller. While we 



as Figure 5(1) shows. The brightest flux density feature 
in the Stokes / spectrum is nicely aligned with the only 
suggestive feature in the Stokes V spectrum. The fits 
agree with what can be seen by eye, and suggest a field 
of 21.9 ± 5.9 mG associated with Gaussian 3. None of 
the other components fitted to the overall profile, listed in 
Table [11] produce a fit regarded as significant. The fea- 
ture in Stokes V being fitted is fairly asymmetric; given 
the large reported field, and relatively narrow line, the 
narrow splitting assumption may be breaking down here, 
and the LCP is amplified more strongly than the RCP. 

IRAS F18588+3517: DG01 discovered this OHM, 
which has multiple narrow components in the 1667 MHz 
line, as well as a weak, but clearly distinguished, 
1665 MH z com ponent. 

Figure 5(m) shows the Stokes / and V features for 
this source, and Table [T^] lists the fits. To fit the 
1667 MHz features, we used six components, includ- 
ing one broad component and three narrow components 
where the emission is strongest. We also fit three compo- 
nents to the 1665 MHz line, which at first glance seems 
a large number given the strength of the 1665 MHz fea- 
tures. Gaussian 2, in particular, is fit to a feature that is 
of smaller amplitude and width than we would ordinar- 
ily fit, but it is motivated by the fact that the velocity of 
the slight peak to which Gaussian 2 is fitted matches ex- 
tremely well with the velocity of Gaussian 5. The widths 
are moderately different, but close enough to plausibly be 
the same feature. The reported magnetic fields are also 
in good agreement, with Gaussian 2 giving 16.3 ± 4.2 
mG, and Gaussian 5 giving 12.8 ± 2.1 mG. Combining 
the two gives 13.9 ±1.9 mG for the region that produces 
the matched 1665/1667 lines. The field fit to Gaussian 
4 is also significant, with a field of 18.3 ± 2.6 mG. A 
large field is also reported for Gaussian 1, but the Stokes 
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TABLE 8 

IRAS F09039+0503 Gaussian Fit Parameters 



S v Au vq By 

Gaussian (mjy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



8.98 ±0.63 1480.2093 ± 0.0013 0.0375 ± 0.0031 37904.2 -16.05 ±2.70 

1 5.25 ± 0.25 1480.4762 ± 0.0097 0.2791 ± 0.0264 37843.3 - 3.39 ± 13.57 

2 6.75 ± 1.73 1480.6593 ± 0.0052 0.0754 ± 0.0129 37801.5 -27.55 ± 5.46 

3 3.69 ±0.57 1480.7525 ±0.0260 0.1400 ± 0.0520 37780.3 - 5.00 ± 14.36 

4 8.66 ±0.36 1481.1259 ±0.0175 0.4110 ± 0.0400 37695.2 -17.81 ± 10.73 

5 4.63 ±0.35 1481.5707 ± 0.0368 0.4333 ± 0.0699 37593.9 -46.32 ± 20.02 

6 6.40 ±0.43 1481.9562 ±0.0029 0.0932 ± 0.0077 37506.1 -12.54 ±5.97 



TABLE 9 

IRAS F10378±1108 Gaussian Fit Parameters 





S 


V 


Au 


Vq 


B \\ 


Gaussian 


(mjy) 


(MHz) 


(MHz) 


(km s- 1 ) 


(mG) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



1.75 ±0.20 1465.5685 ±0.0197 0.4041 ± 0.0575 41277.7 -99.68 ± 43.72 

I 3.56 ±0.37 1466.3600 ± 0.0043 0.0868 ± 0.0110 41093.6 0.56 ± 10.05 

2 2.43 ± 0.54 1466.5915 ± 0.0043 0.0401 ± 0.0106 41039.8 - 3.31 ± 10.24 

3 10.43 ±0.40 1467.1249 ±0.0015 0.0804 ± 0.0039 40915.9 -13.87 ±3.31 

4 17.46 ±0.37 1467.3636 ± 0.0310 3.6798 ± 0.1391 40860.4 -10.08 ± 17.20 

5 4.13 ±0.55 1467.4195 ± 0.0066 0.1265 ± 0.0201 40847.5 20.81 ± 10.92 

6 16.14 ±0.76 1467.5052 ±0.0120 0.7169 ± 0.0195 40827.6 -28.03 ± 6.95 

7 14.99 ±0.74 1467.6153 ±0.0011 0.0730 ± 0.0039 40802.0 4.20 ± 2.21 

8 9.86 ±0.84 1467.7385 ±0.0090 0.2996 ± 0.0237 40773.4 20.95 ± 7.09 

9 2.70 ±0.28 1468.3149 ± 0.0082 0.1569 ± 0.0210 40639.7 9.23 ± 17.79 

10 2.10 ±0.28 1469.5167 ±0.0100 0.1596 ± 0.0255 40361.3 22.77 ±22.98 

II 5.19 ±0.30 1471.2592 ±0.1248 2.9057 ± 0.4441 39958.5 61.95 ± 70.23 



TABLE 10 

IRAS F16255±2801 Gaussian Fit Parameters 



S V Au Vq B|| 

Gaussian (mjy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



8.46 ±2.45 1470.5518 ± 0.0079 0.0528 ± 0.0096 40121.9 - 3.45 ± 3.87 

1 21.22 ±0.82 1470.6097 ± 0.0030 0.0657 ± 0.0115 40108.5 7.01 ± 1.91 

2 16.63 ±1.20 1470.6763 ± 0.0020 0.0482 ± 0.0044 40093.1 9.61 ± 2.02 

3 9.58 ±0.65 1470.7556 ± 0.0023 0.0771 ± 0.0073 40074.8 26.78 ±4.02 

4 6.60 ±0.60 1470.8160 ± 0.0091 0.5459 ± 0.0264 40060.8 64.13 ±14.37 

5 5.99 ±0.52 1470.9155 ±0.0025 0.0824 ± 0.0088 40037.9 10.00 ± 6.07 
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IRAS 18368±3549 Gaussian Fit Parameters 
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V spectrum has noisy features of comparable width, and 
the reported field is not regarded as meaningful. 

IRAS F22134+ 0043: T his O HM was discovered 
with the GBT b y iWillettl (12012ft . T he heliocentric 
rcdshift is 0.212 ([Stanford et all I2000D . but the peak 
1667 MHz emission occurs at a redshift of 0.211. The 
Stokes I spectrum features broad emission mostly be- 
tween 10-20 mjy, and one narrow peak that reaches 
50 mjy. There is possible 1665 MHz emission corre- 
sponding to the 1667 MHz emission redward of the peak, 
but the low amplitude and lack of any 1665 MHz emis- 
sion that corresponds to the 1667 MHz peak suggest that 
the broad feature is bandpass structure. 

Our fit to the emission from this source is shown in 



Figure 5(n) with parameters listed in Table H3l To ade- 
quately capture the broad emission and complex peaked 
structure in the Stokes / spectrum required a total of 
eleven Gaussian components. Fitting the peak was par- 
ticularly difficult. The primary Stokes V feature is an 
asymmetric "S" at the same frequency as the Stokes I 
peak, and producing what was ultimately only a decent 
fit to both spectra at that frequency required three Gaus- 
sian components all centered there, with three very dif- 
ferent widths. These three Gaussians, labeled 5, 6, and 
7, each report very similar magnetic fields: -12.2 ± 4.7, 
-7.1 ± 1.7, and -13.5 ± 2.6. While the field reported with 
Gaussian 5 is not, on its own, considered significant given 
the associated error and width of the line, it is encour- 
aging that all three lines agree so well, given that they 
likely are produced by the same masing clouds. 

In addition to these reported fields, there were also 3a 
fields fitted to Gaussians 2, 8, and 9. As Gaussian 2 is 
wide, it is not considered significant. Gaussians 8 and 
9 are moderately narrow, and the quality of the fits are 
good, with fields of -26.8 ± 5.3 mG and -22.1 ± 4.6 mG, 
respectively. The features being fit are small amplitude 
relative to those associated with the peak emission, but 
are larger than Stokes V structures unassociated with 
Stokes / emission. Unfortunately, shifting the Stokes / 
and V are relative to each other does not provide any 
added support, as the brighter features dominate. So, 
despite the high quality of the fit, we consider the fields 
associated with Gaussians 8 and 9 to be strongly sugges- 
tive, but not confident detections. 

5.2. Linear Polarization 

As in R08, we derived the linear polarized intensity for 
each source from the Stokes Q and U spectra. Unfortu- 
nately we did not find believable linear polarization in 
any of the newly observed OHMs. For this reason, we do 
not provide discussion of the process, instead referring 
the reader to R08. 

6. DISCUSSION 

This work confirms the most basic result of R08, show- 
ing that OH megamasers are viable targets for detection 
of extragalactic Zeeman splitting, and direct measure- 
ment of line-of-sight magnetic fields. Six of the eight 
sources observed by R08 were re-observed here, and the 
Stokes V features and derived magnetic fields were quite 
consistent for five of the six sources. The exception, 
IRAS F10173+0829, featured weak interference in the 
Stokes V spectrum that coincided with the peak of Stokes 



/ emission, at a level barely distinguishable from noise in 
the Stokes / spectrum for the source, and masqueraded 
as Zeeman splitting. 

The success rate in detecting Zeeman splitting 
was significantly lower than in R08. Not counting 
IRAS F10173+0829, four of the eight sources R08 ob- 
served displayed Zeeman splitting, compared to eleven 
of the seventy-one newly observed sources. This is not 
surprising, as R08 selected the highest peak flux density 
OHMs for observations. The peak-to-trough Stokes V 
amplitude of Zeeman splitting features that are less than 
the line width will generally be much smaller than the 
Stokes / feature to which it corresponds. High flux den- 
sity lines are necessary then for detecting weak Zeeman 
splitting. For this reason, all of the newly detected mag- 
netic fields presented here are above the median magnetic 
field amplitude in R08, which was ~3 mG. Of newly de- 
tected fields alone, the median magnitude of magnetic 
fields was ^16 mG, and the median of the sample, count- 
ing re-detections and new detections, is ~ 12 mG. The 
magnetic field fits and errors, as a function of the Stokes 
/ flux density, is shown in Figure [6] 

The detections from R08 are clustered in the lower 
right hand corner, while the new detections mostly fall 
along a trend that reflects that the weaker the flux den- 
sity of lines, the stronger field required for a Zeeman 
splitting detection. The lowest flux density point, at 
3 mjy, corresponds to the lone 1665 MHz Zeeman split- 
ting detection, which nicely agreed with a 1667 MHz 
detection at the same velocity. Many low flux density 
sources are omitted from Figure [H as they did not pro- 
vide meaningful upper limits on the magnitude of mag- 
netic fields. 

It is possible that in addition to missing weak mag- 
netic fields, this survey also misses very strong magnetic 
fields present in OHMs. A few of the detected fields 
were strong enough that the assumption of narrow line 
widths was only marginally applicable. This typically oc- 
curred in sources with very complex emission, in which 
it is likely that circularly polarized emission from phys- 
ically unassociated masing clouds blended together. In 
IRAS 04332+0209, the Stokes / and V spectra suggest 
the possibility that the right and left circularly polarized 
components are completely separated, corresponding to 
a 47 mG magnetic field detection. While confirmation 
of this interpretation is not possible without VLBI ob- 
servations, in a source with more complex emission and 
multiple blended components, it is not even possible to 
identify that such a strong field may exist. The absence 
of fields stronger than ~30 mG inferred using our stan- 
dard analysis procedure may reflect a limitation in our 
approach, rather than strongly disavowing the possibility 
of such strong magnetic fields in OHMs. 

6.1. Comparison with Galactic OH masers 

R08 compared their measured magnetic field strengths 
with results fro m OH masers in the Milky Wa y. Combin- 
ing res ults from lReid fc Silversteinl (|1990ft and lFish et al.l 
(2003), R08 noted that in Galactic OH maser Zeeman 
pairs, the mean was mG, within th e 3.3 mG e r ror in 
the distribution. In a similar survey, iFish et al.1 (|2005ft 
detected a total of 184 Zeeman pairs among all OH 
lines. A co mpanion paper pres enting detailed analysis of 
these data (|Fish &: Reidl 120061) shows their distribution 
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TABLE 12 

IRAS F18588+3517 Gaussian Fit Parameters 



S v Au v Q Bu 

Gaussian (mjy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



O a 1.58 ±0.21 1505.8294 ± 0.0662 0.3804 ± 0.1439 31768.9 6.72 ± 30.63 

1~ 2.32 ±0.48 1506.1660 ±0.0328 0.2691 ± 0.0600 31694.8 57.68 ±17.10 

2 2.98 ±0.65 1506.1751 ±0.0032 0.0318 ± 0.0086 31692.9 16.29 ±4.21 

3 9.53 ±0.20 1507.8172 ±0.0057 0.6455 ± 0.0149 31721.0 -17.45 ± 9.92 

4 10.10 ±0.53 1507.8706 ± 0.0016 0.0463 ± 0.0036 31709.2 18.31 ± 2.60 

5 14.82 ±0.48 1507.9445 ± 0.0012 0.0633 ± 0.0033 31693.0 12.78 ±2.07 

6 6.90 ±0.67 1508.0288 ±0.0013 0.0270 ± 0.0032 31674.5 6.05 ± 2.79 

7 3.50 ±0.41 1508.4066 ± 0.0077 0.1553 ± 0.0207 31591.4 -17.25 ± 13.58 

8 2.99 ±0.19 1508.6885 ± 0.0200 0.3726 ± 0.0489 31529.5 -38.22 ± 24.61 



a Velocities and magnetic fields for Gaussians 0, 1, and 2 correspond to values for the 1665 MHz 
line. 



TABLE 13 

IRAS F22134+0034 Gaussian Fit Parameters 



S v Au vq By 

Gaussian (mjy) (MHz) (MHz) (km s" 1 ) (mG) 
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Fig. 6. — Magnetic field detections and upper limits are shown 
relative to the flux density of the 1667 MHz line. For OHMs with 
detected Zeeman splitting, the flux density used is that of the Gaus- 
sian component with which the Zeeman splitting was associated. 
For upper limits, the plotted flux density is the peak flux density 
of the line. The most confident detected fields are black, moderate 
quality detections are green, and the marginal detections are cyan. 
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Fig. 7. — Three distributions of magnetic field strengths are com- 
pared. The black solid line shows the re sults from the survey of 
Galactic OH mascrs bv IFish et all pOOST ). The dashed green line 
shows the distribution of magnetic fields we found in OHMs. The 
dash-dot cyan line represents the expected distribution we would 
have observed in OHMs if magnetic fields in OHMs had comparable 
strength to magnetic fields in Galactic OH masers. 



of magnetic field magnitudes peaked at 4 mG, and fell 
off quickly thereafter. The two strongest observed fields, 
in W51 e2, had magnitudes of 19.8 mG and 21 mG. The 
results from R08 were broadly similar, as the median 
field detected had a magnitude of ~3 mG, and from that 



R08 concluded that when examined on small scales, star 
formation in (U)LIRGs proceeds similarly to star for- 
mation in the Milky Way. The results from this larger 
survey suggest that strong fields may be more common 
in OHMs than in Galactic OH masers. 
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Direct comparison of the two distributions of detected- 
field s trengths is misleading, however, as the Fish et al. 
(2005) magnetic field limits were far lower than the limits 
of this OHM survey. To better compare the two distribu- 
tions, we simulated an experiment where magnetic field 
strengths were drawn from the observed distribution of 
magnetic fields obs erved in Ga l actic OH masers. This 
was made up of the iFish et all (|2005l ) distribution, plus 
a -40 mG field discovered in W75N (jSlvsh fc Migenesl 
2006). The limits for the simulated survey were drawn 
from our own results. 

The sample of limits was created in the following way. 
For sources with detections, we counted each detected 
field associated with a Gaussian component, and used 
four times the error associated with that component as 
the limit, since our reported upper limits were found by 
taking the reported field strength and adding three times 
the reported error. For sources with no detection, we 
used the upper limit for magnetic fields in the source re- 
ported above. Then, we randomly drew a magnetic field 
strength from the Galactic OH maser magnetic field dis- 
tribution and a limit from our distribution, and recorded 
fields larger than the limits as detections. For a single 
run, this process was repeated n times, where n was the 
number of values in the upper limit distribution. 

We then simulated 10 4 surveys in which the magnetic 
field strengths in the OHM sample were distributed ac- 
cording to the distribution observed in the Milky Way. 
The results of this simulated survey are shown in Fig- 
ure [7J alongside the dis tribution of fields in Milky Way 
OH masers reported in IFish et al.l (|2005l ). and the dis- 
tribution we actually observed. The three distributions 
are clearly distinct. If the magnetic fields in OHMs were 
drawn from the same distribution as those in Galactic 
OH masers, we would have detected many fewer fields 
overall, and, most notably, many fewer strong fields. 
Fields stronger than 10 mG are apparently much more 
common in OHMs than in the Milky Way. 

The results of this simulated survey were consistently 
different than what we observed. The median of the me- 
dian magnitude magnetic field detected in the simulated 
surveys was 5.0 mG, and 95% of the medians fell between 
4.0-5.5 mG. The 12 mG median of OHM magnetic fields 
in the actual survey is more than a factor of two larger 
than the typical median in the simulated samples. Al- 
together, these results argue that the magnetic fields in 
OHMs are stronger than those in Milky Way OH masers. 

6.2. Strong Fields and Time Variability 

R08 reported a factor of —2 variation in one compo- 
nent of the emission from IRAS F12032+1707, and noted 
that the magnetic field associated with that component 
was the strongest observed in their observations. R08 
noted that a connection between strong fields and time 
variability had been ob serve d in Galactic OH m asers, as 
ISlvsh fc Migenesl (I20O6T) and IFish fc Reidl (12001 both re- 
ported 40 mG fields in W75N that are associated with 
maser spots in the process of flaring. 

Our more recent observations of magnetic fields in 
OHMs do not provide any further evidence of this con- 
nection, though the observations do not provide strong 
contrary evidence either. There is weak variability, at 
the 1-10% level, observed in four of the six OHMs re- 
observed in this survey. The strongest variability ob- 



served here occurs in IRAS F12112+0305, in which the 
upper limit on magnetic field magnitude is —3 mG, and 
in III Zw 35, which has magnetic fields with magnitudes 
of —3 mG in the region where variability occurs. In 
IRAS F12032+1707, any further variability since the ob- 
servations of R08 is limited. In comparing other sources 
with strong fields detected in our survey with existing 
spectra in the literature, we encounter no instances sim- 
ilar to IRAS F12032+1707, where a component with a 
strong magnetic field has clearly varied relative to its 
previously published strength. 



6.3. Dynamical importance 

The small number of OHMs with VLBI observations 
show reasonable diversity, with compact, parsec-scale 
emission providing the majority of emission in some 
sources (Arp 220, III Zw 35, IRAS F17208-0014), while 
compact emission is weak or absent in others (Mrk 231, 
IRAS F14070+0525), and in one case, emission was 
compact on the ~100 pc resolution of the observations 
(IRAS F12032+1707). Of the strong Zeeman splitting 
signals observed in Arp 220 and III Zw 35, most could 
be associated with these compact, parsec-scale features. 
If all strong magnetic fields detections presented here, 
including those in sources without VLBI observations, 
occur in parsec-scale regions, it suggests that magnetic 
fields may be dynamically important in masing clouds in 
OHMs. 

In the m o deling of OHM emission in III Zw 35 by 
iParra et"al~l (|2005f) . they assumed a maximum density 
of uh 2 = 10 5 cm~ 3 in their OHM clouds, a cloud size 
of roughly R = 1 pc, and noted that typical internal 
velocity dispersions are AV — 20 km s _1 . We take 
these values as being reasonably representative of OHMs 
at large, and explore the implications of the observed 
magnetic field strengths. While this is surely imperfect, 
the small sample of VLBI measurements indicates that 
III Zw 35 is no t unusual in its small s cale properties. 

Following the lStahler fc Pallal (|2005l ) application of the 
virial theorem to a spherical cloud, and assuming its mass 
is dominated by molecular hydrogen, the ratio of the 
magnetic energy density, M, to the self gravitational en- 
ergy, W, is 



M 

|W| 



0.5 



D 



3 mG 



R 
1 pc 



10 5 



,-3 



(4) 



R08 found —3 mG line-of-sight magnetic fields in 
III Zw 35. For a variety of plausible magnetic field proba- 
bility distribution functions, the mean of the line-of-sight 
magnetic fields is ro ughly half the mean of th e true mag- 
netic field strength (jHeiles fc Crutcherl 120051) . This sug- 
gests that the typical magnetic fields observed in OHMs 
a re dynamical l y imp ortant, and perhaps even dominant. 

IParra et al.l (|2005l ) note, however, that the clouds in 
III Zw 35 have typical velocity dispersions of —20 km s~ 
and are not likely gravitationally confined, and then go 
on to suggest that magnetic fields could play a role in 
magnetically confining the clouds. The ratio of the mag- 
netic energy density and the turbulent energy density, T, 
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are equal for 



M 

\r\ 



0.2 



B 



3 mG 



AV 



20 kms" 



10 5 cm -3 



(5) 

This expression shows that for total magnetic field 
strengths a factor of two stronger than the observed linc- 
of-sight fields, an d values for t he ve locity dispersion and 
cloud size used in lParra et al.l (|2005f ) . the turbulent pres- 
sure and magnetic pressure are of roughly equal strength. 
Altogether, these results argue that magnetic fields are 
important in determining the structure of OH masing 
clouds in at least some fraction of OHMs. 

Densities in Galactic OH ma ser regions are of order 
10 6 -10 7 cm" 3 (jReid et alJll987f >. as measured by NH 3 . 
If extragalactic OH megamasers mase at similar densi- 
ties, the weaker (~3 mG) magnetic fields detected in 
this survey would not be dynamically important, unless 
clouds are a factor of a few smaller than 1 pc. However, 
for sources with ~20-30 mG fields, magnetic fields would 
still be dynamically relevant in a ~1 pc cloud, even at 
densities of 10 6 cm" 3 . More VLBI observations of OHMs 
with magnetic fields, particularly the strongest magnetic 
fields, will be necessary to developing a better under- 
standing of the role magnetic fields play in the central 
regions of OHM hosts. 

7. CONCLUSIONS 

We performed Full-Stokes observations of seventy- 
seven OH megamasers, representing the entire known 
sample of Arecibo accessible OHMs known at the time 
of observations. Six of these so urces had recent l y bee n 
observed in Full-Stokes mode by Robishaw ct al. (2008), 
and four of those six sources had features in their 
Stokes V spectra consistent with Zeeman splitting. We 
confirm the Stokes V detections in three of the four 
sources, the non-detections in the two sources without 
Stokes V features, and find that the apparent Zeeman 
splitting in one of the four sources was actually inter- 
ference. Of the remaining seventy-one sources without 
previous Stokes V observations, eleven have features con- 
sistent with Zeeman splitting in the Stokes V spectra. 



For all sources with Stokes V features, we derived 
magnetic field strengths that could produce the observed 
Zeeman splitting signal. The median magnetic field 
magnitude associated with masing components in the 
fourteen OHMs with Zeeman detections is 12 mG. 
We show that the magnetic field strengths observed 
in OHM masing clouds are roughly a factor of 2-3 
larger than those observed in Galactic OH masers. 
The strongest derived magnetic field in which we are 
confident has a magnitude of 27.6 mG. One source, 
IRAS F04332+0209, appears to have Zeeman splitting 
that exceeds the linewidth. This interpretation cannot 
be confirmed without VLBI observations, but if it is 
correct, the required magnetic field strength would be 
~47 mG. In two other sources, IRAS 02524+2046 and 
IRAS FI8588+35I7, we observed for the first time 
Zeeman splitting in the 1665 MHz OH line of an OHM. 
In each case, the reported magnetic field strength is 
consistent with that derived for the 1667 MHz line at 
the same velocity. For reasonable assumptions about 
conditions in OHMs, the magnetic field measurements 
presented here argue for magnetic fields playing a 
dynamically important role in the masing clouds. 
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